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ABSTRACT 

Aims. We search for young stellar objects (YSOs) in the LI 64 IN cluster and characterize the star formation activity through determination of 
the age distribution, mass function, spatial distribution, and the star formation history. 

Methods. Multi-wavelength broad band photometry both from space and the ground are used to look for IR excess in order to separate field 
stars from YSOs and to sample the spectral energy distributions. Space-based observations were obtained using the ISO satellite (ISOCAM) 
in two filters, centred at 6.7 and 14.3 /jm, and Spitzer (IRAC) at 3.6, 4.5, 5.8, and 8.0/jm. Our ground-based observations were made with the 
Nordic Optical Telescope (NOT) using ALFOSC (/ band), NOTCam (/, K s and 2.12pm H 2 ), and SIRCA (U). More than 50 of the brightest 
/-band sources were then studied with follow-up optical spectroscopy (5780-8340 A) to check for signs of accretion (Ho- in emission) and 
youth (Li I /16707 in absorption) and to determine their effective temperatures. By comparing theoretical evolution tracks with our YSO sample 
in the H-R diagram, we calculated an age, luminosity, and mass distribution. 

Results. We detect a total of 216 (Spitzer or / band) sources in L1641N, 89 of which are YSO candidates. Most of the spectra are of M-type 
with Ha strongly in emission, and many have Li 6707 in absorption. The four brightest / band sources (F and G stars) are suggested as 
foreground stars, and the L1641N IRAS source is shown to be the combined flux of at least four sources. We find that the interstellar extinction 
is well-fit in the optical and near-IR by a power law with an exponent of 1.58, although in the mid-IR the Spitzer observations show a higher 
extinction than expected from theory. The median age of the YSO sample is — 1 Myr and the resulting MF has a flat distribution for low masses 
down to the completeness limit. There is evidence of a constant star formation rate of one star in 3.7 x 10 4 yr during the past few Myr. We find 
11 sources older than 10 Myr and a spatial separation between younger and older YSOs, suggesting that many of the older stars formed in 
L1641N could have left the cluster, giving the appearance of an increased star formation rate with time. 
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1. Introduction 

Lynds 1641 (L1641) is a dark molecular cloud, located in the 
southern part of our nearest giant molecular cloud Orion A 
(d~450pc). It is a site of active star formation with a lot of 
very young stars. L1641 is therefore ideal for studying both the 
stars themselves and the related Herbig-Haro (HH) outflows. 
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Attention was drawn to the northern part of the cloud, 
simply called L1641N, in the mid 80s with the discovery 
of a molecular outflow (Fukui et al. 1 19861 1 at the position 
of the mid-IR bright source IRAS 05338-0624. Many stud- 
ies have investigated the nature of this source, Fukui et al. 
(1988 ) found well-separated lobes to the north (blue-shifted) 
and south (red-shifted) of the IRAS source using CO observa- 
tions of higher resolution. Further observations have been made 
through molecular line studies (e.g. Chen et al. 1996 Sakamoto 
et al. 119971 Stanke & Williams |2007>, Optical and near-IR sur- 
veys (e.g. Strom et al. [19891 Hodapp & Deane [T9931 Chen et 
al. [T99"3l and in the mid-IR (Ali et al. l2004t . 

L1641N is presently the most active site of low-mass star 
formation in the L1641 molecular cloud and has one of the 
very highest concentrations of HH objects known anywhere in 
the sky (Reipurth el al. 1998). HH shocks in these flows have 
been observed in both the near-IR 2.12 fim line of H2 (Stanke et 
al. [1998 ) and in the usual optical shock lines Ha and the [SII] 
doublet (Reipurth et al. [T998l and Mader et al. fl999b . 
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In this article we make the most detailed mid-IR study yet 
of L1641N using both our ISO observations and the much 
more sensitive and higher resolution Spitzer data to survey 
the region for YSO candidates and determine the mass func- 
tion in L1641N. Even though the Spitzer observations cover 
a much larger region, we have restricted ourselves to the re- 
gion we selected for our ISOCAM observations in order to get 
a more complete wavelength coverage and because L1641N is 
the most dense, active part of L1641. We have also made ad- 
ditional near-IR and optical ground-based imaging and taken 
optical spectra of many of the objects, to be able to confirm the 
YSO status and to estimate their effective temperatures. Finally, 
as part of this study, we have also made deep narrow band ob- 
servations using a 2. 12 /mi H2S(1) filter to study the numerous 

sources in the region (Galfalk & Olofsson 120071) . 
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Fig. 1. System transmission curves of IRAC channels 1-4 
(Spitzer). The four major bands of PAH emission are marked 
by dotted lines. IRAC channel 4 (centred at 8.0 /mi) contains 
two strong bands while images in channel 2 have very smooth 
backgrounds as they are essentially free from PAH emission. 



2. Observations and reductions 

2.1. Space based 

Mid-IR observations were carried out using six filters centred 
in the wavelength range 3.6 - 14.3 /mi using ISOCAM on- 
board the ISO satellite and IRAC onboard the Spitzer Space 
Telescope. 

2.1.1. ISO satellite 

We have used the ISOCAM instrument onboard the ISO satel- 
lite (60 cm primary mirror) to observe L1641N at mid-IR wave- 
lengths as part of a mid-IR star formation survey (princi- 
pal investigator L. Nordh). Our observations were made us- 
ing two broad-band filters at 6.7 and 14.3//m. Centred on the 
IRAS source 05338-0624 they cover an area of 835 x T.60 
~ 0.0176 sq.deg. (63.5 sq.arcmin). 

All imaging was done at a PFOV of 3" and the integra- 
tion time of each transmitted 32 x 32 pixel frame was 2.1 s. 
For both filters the full 10x9 mosaic of 90 images uses 1550 
frames, thus at least 17 frames are combined at each position, 
forming images with a temporal history which is used to e.g. 
remove cosmic rays and source transients (ghost sources) from 
previous mosaic positions. There is also a spatial redundancy 
in both equatorial directions due to image overlaps, which is 
necessary to avoid having a source imaged on the same array 
pixels, especially on the unfortunate dead column (No 24). This 
of course also means that we have an increased total exposure 
time and temporal history in overlapping regions. 

For the data reduction we used the CIA V4.0 package 
(Ott et al. [19971 Delaney et al. EOOOi l and the SLICE package 
(Simple & Light ISOCAM Calibration Environment) accessed 
from inside CIA. The CIA reduction steps consisted of (in or- 
der): Extracting useful observations, dark correction (Vilspa 
dark model, Biviano et al. 2000), glitch removal (multireso- 
lution median transform, Starck et al. 1 19961 1. short transient 
correction (Fouks-Schubert model, Coulais & Abergel 2000), 
flatfielding (constant median flatfield from observations). We 
also used a set of our own programs for source detection and 
photometry. 



As a last step, all frames were projected into mosaics in 
both filters, followed by point source detection and aperture 
photometry. All point sources were traced back to their corre- 
sponding original frames for temporal and spatial verification 
(in order to exclude remaining artefacts). For the photometry, 
the point spread function was used to correct for the flux out- 
side of each aperture. In total 41 sources were detected, with 
photometry possible for all sources at 6.7 /an and for 31 of 
these at 14.3 /mi. 

Even though the integration time was 2.1 s for each frame 
and there are about 17 frames for each image, the total expo- 
sure time for a given pixel in the mosaic varies between 36 s 
and 143 s due to the varying number of overlaps. Thus, it is ex- 
pected that the faintest sources will generally be detected away 
from the mosaic edges. From mean temporal and sky noise 
calculations of all detected ISOCAM sources in L1641N, the 
photometric uncertainties are estimated to be (lcr) 2.9 mJy at 
6.7 /mi and 5.4 mJy at 14.3 /mi which also approximately rep- 
resent the detection limits. However, such faint sources are only 
detected in low nebulosity regions that are free of artefacts. On 
the other hand the surveyed region is not completely mapped 
down to the lcr level, since it is possible that even sources 
brighter than 3cr (8.7 mJy and 16.2 mJy respectively) may have 
been unnoticed due to varying nebulosity, glitches, memory ef- 
fects, uncovered dead columns and source confusion close to 
very bright sources. 

In order to convert mJy fluxes into 6.7 //m and 14.3/mi 
magnitudes, we have used the zero-magnitude flux densities 
given in Blommaert et al. (2003) which gives the following re- 
lations: 

m 6 . 7 = 12.39 - 2.51og 10 F 6 . 7 
m 14 . 3 = 10.74 - 2.5 logi F14.3 

2.1 .2. Spitzer Space Telescope 

Spitzer is an IR space telescope that carries a 85 centimeter 
cryogenic telescope and three cryogenically cooled science in- 
struments, one of these is the Infrared Array Camera (IRAC) 
that provides simultaneous 5'.2 x5'2 images in four channels, 
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Table 1. Spitzer archive data used. 



Key 


Type 


Released 


Scheduled 


4101888 


IRAC Mapping 


2005-07-27 


2004-02-16 


4102912 


IRAC Mapping 


2005-07-27 


2004-02-17 


4101376 


IRAC Mapping 


2005-07-27 


2004-02-18 


4102400 


IRAC Mapping 


2005-07-27 


2004-02-18 


4101632 


IRAC Mapping 


2005-10-27 


2004-10-08 


4102656 


IRAC Mapping 


2005-10-27 


2004-10-08 


4102144 


IRAC Mapping 


2005-11-17 


2004-10-27 


4103168 


IRAC Mapping 


2005-11-17 


2004-10-27 



Table 2. Spitzer zero-magnitude flux densities. 



Channel 


F(0) 


[Jy] 


Channel 1 (3.6 /im) 


280.9 


±4.10 


Channel 2 (4.5 jum) 


179.7 


+ 2.60 


Channel 3 (5.8 yum) 


115.0 


± 1.70 


Channel 4 (8.0 fim) 


64.13 


±0.94 



centred at 3.6, 4.5, 5.8 and 8.0 fim (Fig. [T}. Each channel is 
equipped with a 256 x 256 pixel detector array with a pixel 
size of about l'/2 x 1'.'2. 

The Spitzer data used in this paper (see Table [TJi was ob- 
tained from the Spitzer Science Archive using the Leopard soft- 
ware. All the data is part of the program "An IRAC Survey of 
the L1630 and L1641 (Orion) Molecular Clouds" (Prog.ID 43) 
with G. Fazio being the P.I. All data had been reduced to the 
Post-Basic Calibrated Data (pbcd) level, however since there 
was still a lot of artefacts, background variations and varying 
orientation in the mosaics we used our own in-house routines 
to reduce the data further. This was possible since we had eight 
overlapping mosaics to merge in each channel, including cos- 
mic ray removal, de-striping (most likely pick-up noise), re- 
moving the large background variations between regions in the 
same mosaic and also between all mosaics, marking bad pix- 
els not to be used further, removing ghost-effects from bright 
sources, de-rotation, shifting, adding and making a compos- 
ite colour image. This required some work, but since we only 
needed to reduce the L1641N region (and these mosaics are 
huge) we could concentrate on just that part - which because 
of all mosaics having different rotations meant first finding 
L1641N in each one and then cropping them to manageable 
sizes. 

Table [2] gives the zero-magnitude flux densities, F(Q), that 
have been used for the flux-magnitude conversion in the four 
channels (also see e.g. Reach et al. 120051) . The conversion for a 
measured flux Ft in channel i is then m, = 2.5 log\o(F(0)/Fj) 

2.2. Ground based 

All photometry and spectroscopy were carried out at the 2.56 m 
Nordic Optical Telescope (NOT) located at 2382 m (7815 feet) 
above sea level on the island of La Palma, Canary Islands. 
We used three different cameras: ALFOSC (Optical), NOTCam 
(Near-IR), SIRCA (Near and mid-IR). For details, see the sub- 
sections below. 



2.2.1 . Optical photometry (/ band) 

The / band observations were obtained on Dec. 04 2003 us- 
ing the ALFOSC (Andalucia Faint Object Spectrograph and 
Camera) with an average seeing of about l'.'l (in the / band). 
This instrument has a 2048 x 2048 CCD and at a PFOV of 
07188/pixel it has a FOV of about 6'4 x 6!4. A 5 x 5 mosaic 
was made with individual exposure times of 60s using step 
sizes of 23" and 30" in RA and Dec respectively, meaning a 
total exposure time of 25 minutes throughout most of the mo- 
saic. Aperture photometry was then carried out on all visible 
sources using our own ALFOSC IDL-scripts. 

2.2.2. Optical spectroscopy (5780-8340 A) 

ALFOSC long-slit spectra were obtained during six nights at 
the NOT, Dec. 02-04 2003 and Jan. 08-10 2005, using grism 
#8 which yields spectra in the wavelength range 5780-8340 A. 
We used a fairly wide slit (1"2) in order to minimize slit losses 
(the seeing was on average ~l'/0). 

The selected wavelength range covers some very interest- 
ing lines for low-mass star formation, as they are signatures of 
youth (Li 6707 in absorption), accretion (Ho- in emission) and 
a number of TiO and VO bands that can be effectively used in 
the spectral classification of late-type stars. Since these sources 
are relatively faint in the optical we had to use total exposure 
times ranging from 40-120 min for many sources, but when- 
ever possible we tried to put 2-3 stars on the same slit. The 
total exposure time used for the spectroscopy is about 35 hours. 

Standard subtraction of bias and dark current was carried 
out, however, the spectroscopic flatfielding was complicated by 
strong fringing. Great care had to be taken in order to correctly 
remove this fringing, which otherwise becomes a problem in 
the red part of the spectrum. Halogen flats were taken at each 
source position directly after each observation since the pat- 
tern varies with position and time. Also, a constant was fit and 
added prior to flatfielding in order to fine-tune the pattern re- 
moval because of background variations. For full-size frames 
this constant was found from the overscan region, and iterated 
in sub-frame exposures. 

For wavelength calibration we used 29 Ne lines in long and 
short exposures of a Ne calibration lamp to be able to measure 
both bright and faint lines. A polynomial fit with all these Ne 
lines included was made at each spatial position along the slit 
(meaning 2048 solutions). These fits were also used to unbend 
the spectra in the wavelength direction. Distortion correction 
was made in the spatial direction using calibration stars at dif- 
ferent positions on the slit. For each spectrum a constant shift 
was also applied to correct for small shifts of the whole wave- 
length scale with different alt-az positions of the telescope (for 
this we used several skylines). 

2.2.3. Near-IR photometry (/, K s , V and 2.12/im H 2 ) 

The J band images were observed with NOTCam (near-IR 
camera) on the night of Dec. 10 2003 using 70 exposures of 
48 s with random small-step dithering in between. These ob- 
servations cover the central part of LI 64 IN as illustrated by the 
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Table 3. L' band photometry of the central region. 
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Fig. 2. SIRCA L' band image of central L1641N. The field 
shown has a size of 73" x 47". 



green square in Fig. [3] (NOTCam has a square FOV of 4' and 
a pixel size of 07235) and was made as a complement to the 
2MASS / photometry to measure fainter sources. There were 
thin cirrus covering the sky that night so we have used relative 
photometry with 2MASS sources for calibration. 

We have selected two regions for K$ band imaging (48 and 
24 minutes total exposure times, respectively). These obser- 
vations were made together with deep 2. 12 /mi H2 S(l) ob- 
servations (175 and 100 minutes) as part of a related project 
to measure the proper motions of H2 objects in L1641N 
(Galfalk & Olofsson 2007). In the present paper we focus 
on the Ks photometry, for which we detect point-sources 
down to Ks ~ 19.5 mag. The observations are also used to 
detect wide double sources and for aperture photometry of 
these. The observations were made on two photometric nights, 
Dec. 13-15 2005, with an average seeing of 0775 (0760- 
0785) using NOTCam with the newly installed science-grade 
array. NOTCam is an HAWAII 1024 x 1024 x 18.5 jum pixels 
HgCdTe array with a field-of-view of 4!0 x 4!0. Both regions 
include the centre of L1641N, which thus has an increased 
total exposure time. The first field is centred on a position 
(05 ft 36 ffl 24.07*, -06°23'01'.'9, Epoch 2000) close to the bright- 
est mid-IR source in L1641N (No. 172 in this paper). The 
second field, centred on (05''36 m 13.94 s , -06°20'52'.'5, Epoch 
2000) images the NW and central part of L1641N. 

Besides the usual reduction steps of near-IR imaging, we 
have used our NOTCam model (Galfalk 2005) to correct for 
image distortion and some other in-house routines (written in 
IDL) to find and remove bad pixels, shift-add images and to 
remove all the dark stripes that results from lowered sensitivity 
after a bright source has been read out of the detector. 

To probe the central ~ 1' of LI 64 IN (Fig. |2]i we have used 
our own instrument SIRCA (Stockholm IR Camera) for obser- 
vations in the L' band (3.8/mi). Observations were carried out 
20-21 December 2002. SIRCA uses a 320x256 pixel InSb ar- 
ray, yielding a FOV of 70" x56" when mounted on the NOT. In 
these observations we used an individual exposure time of 200 
ms, 5 readouts per chopper position, 5 chopper cycles and 54 
ABBA nodding cycles of the telescope. This gives a total expo- 
sure time of 0.2x5x2x5x4x54 = 2160 seconds (36 minutes). 



INO. 


L (mag) 


Comment 






110 


11.27 ±00.05 


Double source 






117 


13.50 ± 00.43 








123 


12.34 ±00.13 


Deeply embedded (Ks 


-U 


= 5.05) 


124 


09.56 ± 00.02 


Deeply embedded (K s 


- L' 


= 6.84) 


126 


14.8: 








134 


13.45 ± 00.40 








138 a 


09.66 ± 00.02 








138 b 


10.44 ± 00.03 









The weather was not fully photometric that night so we had 
to use relative photometry from a previous test night (which of 
course had a much shorter total exposure time) to be on the safe 
side. Therefore the uncertainties may seem larger than the ex- 
posure time would suggest. Flatfielding was done differentially 
by chopping between the cold interior of the camera and the 
much hotter mirror cover of the telescope. 

2.2.4. Additional observations 

We have used the 2MASS archive (Cutri et al. 12003-2006] ) for 
regions of LI 64 IN where we lacked J or Ks photometry. Given 
that this is an all-sky survey, it cannot have the high sensitiv- 
ity and resolution obtainable from target specific observations. 
It still does a fairly good job filling in photometry of sources 
outside of our selected J and K s fields. 

3. Results and discussion 

3.1 . General impression and mosaics 

The final reduced Spitzer mosaic (using our own routines) is 
presented in Fig. [3] as a colour image using three of the four 
IRAC channels (blue = 3.6 /mi, green = 4.5 /mi and red = 
8.0 /mi). One outstanding feature in this image is the very bright 
extended emission at 8.0/mi (also present in the 5.8 /mi chan- 
nel but not as bright). The reason for this emission can be seen 
in Figure[T| it is caused by Polycyclic Aromatic Hydrocarbons 
(PAHs), molecules built up of benzene rings which emit most 
strongly in IRAC 4, followed by IRAC 3, IRAC 1 and finally 
IRAC 2 which is essentially PAH free. 

Since the PAH emission in LI 64 IN is very bright and seen 
across the whole mosaic (see Fig. [3]) the UV radiation causing 
the excitation must be illuminating the cloud uniformly. We 
suggest that the illuminating source is l Ononis (HD 37043), lo- 
cated about 30' to the north. This star has a 09 III spectral type 
and a parallax of 2.46 ± 0.77 mas according to the Hipparcos 
catalogue, corresponding to a distance of about 400 pc (310- 
592 pc). This fits in very well with the assumed distance to 
L1641N of about 450 pc and puts the suggested star at a similar 
or somewhat closer distance than L1641N. Their separation in 
the plane of the sky, if equal distances are assumed, is about 
3.5pc. 
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Fig. 3. Spitzer composite of 3.6/im (blue), 4.5 fim (green) and 8.0/vm (red). All our sources are marked by circles. Sources for 
which we have spectroscopy have thick, dashed circles. The surveyed region has a size of 8f45 x 7!38 (1.11 X 0.97 pc at a distance 
of 450 pc), although the field size in this Figure has been made larger to give a better overview of the region. The green square 
outlines our J band observations. 



Our two K s mosaics (and shocked 2.12yum H2) are pre- 
sented in Figures 22 and 23 with all K s counterparts to Spitzer 
sources marked by circles. There is a large amount of shocked 
H2 in these mosaics, much of which can also be seen in the 
IRAC images (especially in the 4.5 /itn channel). For a further 



discussion on the H2 flows see Galfalk & Olofsson (2007) and 
Stanke et al. ( fT998l 120001 . 
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Fig. 4. Comparison of ISO and Spitzer photometry. The 
ISOCAM 6.7/mi channel is compared to the mean flux of 
the Spitzer (IRAQ 5.8 and 8.0/mi channels. The upper panel 
shows our ISOCam photometry and the lower panel that of Ali 
et al. J2004t . 

3.2. Photometry and source catalogue 

We have carried out photometry in all filters of the ground- 
based (/, J, K s , U) and space-based (3.6, 4.5, 5.8, 6.7, 8.0 and 
14.3 jt/m) observations within the area covered by the ISO ob- 
servations. The full source list is basically the full Spitzer sam- 
ple in L1641N (any sources seen with IRAC are included in 
the list), with (a few) additional sources added from the / band 
mosaic. 

In total we include 216 sources in our L1641N catalogue. 
In order to keep the Table width down we have divided our 
photometry into two Tables. Table|9]contains space-based pho- 
tometry from the ISO and Spitzer telescopes and Table 10 the 



ground-based photometry. We have also observed a small num- 
ber of these sources in the L' band (see Table BJ. Source num- 
bering is by a common cloud number, i.e. the same for all filters 
and increasing from west to east. 

Spitzer sources that appear double in our ground based ob- 



servations (cf Sect. 3.7.3 1, will in the following be called dou- 
ble sources, including both gravitationally bound and optical 
binary sources. In all Tables we use bold source numbers for 
sources where we have found intrinsic IR excess. 

One way of checking the ISO and Spitzer photometry, al- 
though observed using different filters, is to plot the ISOCAM 
LW2 channel (6.7 //m) versus the mean of IRAC channels 3 
(5.8 //m) and 4 (8.0/mi). Points should then roughly follow a 
line with fluxes being equal on both axes. This is shown in 
Fig. fusing both our ISOCAM photometry and that of Ali et 



al. ((2004). Our ISOCAM photometry agrees very well with the 
more accurate Spitzer photometry, whereas the ISOCAM pho- 
tometry of Ali et al. seems to systematically underestimate the 
fluxes by a large factor (even more so at 14. 3 /mi) and as an 
average their flux densities needs to be multiplied by a factor 
of 1.41 (6.7 /mi) and 1.88 (14.3/mi) respectively to agree with 
the Spitzer data. 

Our ISOCAM circular aperture photometry uses well tested 
routines which have been used in several previous papers (Kaas 
et al. [T999l Olofsson et al. [19991 Persi et al. l2000l Bontemps et 
al. EOOn Kaas et al.|2004]and Galfalk et al. 120041 our Spitzer 
photometry uses another set of routines that we wrote espe- 
cially for the Spitzer data. That fact that they agree very well 
with each other, although using different routines and calibra- 
tions (not done by the same people), indicates that both our 
ISO and Spitzer photometry are reliable. In total we detect 211 
sources in the Spitzer mosaics and 43 of these in the ISOCAM 
images (Ali et al. 2004 detected 34 sources). 

3.3. The central region and IRAS 05338-0624 

The optically invisible source IRAS 05338-0624 is one of the 
most luminous IRAS sources detected in the region and in fact 
the only IRAS source in L1641N. It has been suggested that 
its counterpart is N15 or Nl of Chen et al. ( 1993), correspond- 
ing to L1641N-124 and 116 in our numbering scheme. Given 
the uncertainty ellipse of the IRAS measurements and the new 
sources we see in our Ks and Spitzer observations, it can be 
seen in Fig.|5]that the situation is probably much more compli- 
cated than just one source being its counterpart. 

Within the IRAS beam we detect four sources: L1641N- 
115, 116, 123 and 124, all fairly bright in the mid-IR, so we 
conclude that the IRAS source is in fact the added flux of a 
number of sources (with No. 124 clearly being the brightest in 
the mid-IR). There is also some mid-IR reflection nebulosity in 
the central region that can be part of this flux, as can be seen in 
our L' image (Fig. [2}. 

Chen et al. ( 1 19951 1 observed a dust peak in the 2 mm con- 
tinuum (3"beam) close to the IRAS position and suggested that 
this could be a circumstellar dust disk. This scenario fits very 
well with our observations, as a jet with a chain of H2 knots 
seems to emanate from this position. Our L' image only shows 
extended nebulosity at L1641N-116 (Nl), suggesting that only 
reflection nebulosity is seen in the 3.6 /mi Spitzer data. This 
could also be the case for the 4.5 /mi channel (possibly also 
some H2 emission) and the M band observations of Chen et 
al. (Q993) and at even longer wavelengths given the ISO and 
Spitzer fluxes (ISO has a larger beam and thus measures more 
extended flux). 

There is one more bright mid-IR source in the vicinity of 
the central region, and that is L1641N-172. However, given that 
it is located more than an arc minute (86") from IRAS 05338- 
0624, and the fact that it has brightened considerably the 
last decade (Galfalk & Olofsson 2007) makes this a very un- 
likely counterpart. Since only one IRAS source was detected in 
L1641N, L1641N-172 was probably much fainter at the time 
of the IRAS observations and thus not detected. 
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Fig. 5. The central region of L1641N in a composite using / (blue colours), 2.12 //111H2 S(l) (green) and K$ (red) filters. The big 
ellipse marks IRAS beam centred on IRAS 05338-0624, circles denote our sources and plus signs are corresponding sources of 
Chen et al. d 1993b . The square is centred on the 2 mm dust peak (see text). The field shown has a size of 83'.'0 x 8075 for both 
images. 



The source density is clearly much higher at the centre of 
L1641N than in the outer parts, and judging from the H2 flows 



(Figures 22 and 23 1 there should be several very young (Class 
0) YSOs in that region. These are however much too deeply 
embedded to be included in our MF calculations (which needs 
detections in the / and J bands). For further studies of the cen- 
tral region we refer to Chen et al. (1996) and Stanke & Williams 
(20073. 



3.4. Identifying stars with intrinsic IR excesses 

3.4.1 . Intrinsic colours for normal stars 

The main problem of using available empirical colour infor- 
mation for normal stars is the differences between the detailed 
response functions, which may result in relatively large er- 
rors, in particular for late M stars. For this reason we use 
theoretical colours calculated by convolution of the filter re- 
sponse functions and synthetic spectra. We have chosen the 
BaSeL 2.2 spectral library with corrected SEDs (Lejeune et 
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Fig. 6. Stellar models (BaSeL v2.2) for logg = 4.5 and T eff = 2 000 - 10000 K. Overplotted are the ground based filter trans- 
mission curves (7, J and K$ ) used in our observations. 




4 6 8 10 

Wavelength [urn] 

Fig. 7. Stellar models (BaSeL v2.2) for logg = 4.5 and T eff = 2000 - 10000 K. Overplotted are the IRAC filter transmission 
curves. 
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Table 4. Calculated colours of normal stars. The model atmospheres (BaSeL v2.2) and transmission curves used in the calcula- 
tions are shown in Figures[6]and[7] 



T 

T eff 


1-J 


1-K S Kg -3.6 Kg -4.5 
Main sequence, log g = 


1/ CO 


on 


T T 
1-J 


T L' 


Supergiants. 


JVS"4.J 
log g = 


IT Z Q 

Ks-O.5 
U.U 




2400 


2.502 


1.060 


0.646 


0.589 


0.361 


0.502 


3.794 


0.741 


0.003 


0.568 


0.704 


0.838 


2600 


2.244 


0.927 


0.531 


0.497 


0.363 


0.501 


3.136 


1.194 


0.397 


0.775 


0.905 


1.035 


2800 


2.076 


0.846 


0.455 


0.441 


0.377 


0.513 


3.054 


1.270 


0.387 


0.689 


0.807 


0.926 


3000 


1.879 


0.810 


0.371 


0.369 


0.338 


0.466 


2.969 


1.308 


0.329 


0.528 


0.638 


0.745 


3200 


1.647 


0.804 


0.290 


0.263 


0.236 


0.358 


2.524 


1.282 


0.253 


0.320 


0.421 


0.519 


3400 


1.382 


0.822 


0.212 


0.144 


0.163 


0.278 


1.869 


1.188 


0.225 


0.217 


0.304 


0.422 


3600 


1.133 


0.814 


0.173 


0.129 


0.222 


0.325 


1.395 


1.092 


0.186 


0.005 


0.069 


0.246 


3800 


1.031 


0.800 


0.156 


0.103 


0.183 


0.276 


1.146 


0.971 


0.158 


-0.022 


0.032 


0.201 


4000 


0.861 


0.776 


0.151 


0.078 


0.145 


0.229 


1.053 


0.841 


0.139 


-0.033 


0.011 


0.173 


4200 


0.877 


0.717 


0.124 


0.037 


0.093 


0.169 


0.926 


0.745 


0.131 


-0.028 


0.008 


0.161 


4400 


0.871 


0.670 


0.122 


0.032 


0.079 


0.148 


0.858 


0.658 


0.124 


-0.022 


0.009 


0.149 


4600 


0.840 


0.623 


0.118 


0.032 


0.072 


0.135 


0.785 


0.587 


0.119 


-0.012 


0.016 


0.141 


4800 


0.781 


0.572 


0.108 


0.028 


0.063 


0.120 


0.706 


0.531 


0.115 


0.004 


0.031 


0.137 


5000 


0.693 


0.523 


0.099 


0.027 


0.056 


0.107 


0.648 


0.485 


0.112 


0.023 


0.050 


0.135 


5200 


0.627 


0.490 


0.092 


0.029 


0.053 


0.097 


0.601 


0.427 


0.130 


0.068 


0.094 


0.155 


5400 


0.560 


0.443 


0.086 


0.032 


0.052 


0.088 


0.537 


0.378 


0.131 


0.095 


0.119 


0.160 


5600 


0.504 


0.392 


0.080 


0.036 


0.051 


0.080 


0.483 


0.330 


0.125 


0.110 


0.129 


0.157 


5800 


0.466 


0.346 


0.074 


0.040 


0.050 


0.071 


0.449 


0.286 


0.117 


0.116 


0.131 


0.150 


6000 


0.450 


0.316 


0.070 


0.044 


0.050 


0.064 


0.433 


0.257 


0.113 


0.122 


0.133 


0.147 


6200 


0.398 


0.277 


0.062 


0.045 


0.048 


0.056 


0.421 


0.232 


0.110 


0.126 


0.135 


0.144 


6400 


0.347 


0.236 


0.055 


0.044 


0.044 


0.049 


0.351 


0.187 


0.102 


0.122 


0.132 


0.143 


6600 


0.295 


0.197 


0.048 


0.043 


0.040 


0.042 


0.279 


0.144 


0.096 


0.118 


0.129 


0.142 


6800 


0.252 


0.163 


0.041 


0.039 


0.034 


0.034 


0.230 


0.112 


0.091 


0.115 


0.126 


0.139 


7000 


0.230 


0.137 


0.034 


0.032 


0.027 


0.026 


0.203 


0.090 


0.084 


0.108 


0.119 


0.132 


7200 


0.183 


0.127 


0.027 


0.027 


0.022 


0.021 


0.148 


0.082 


0.077 


0.101 


0.112 


0.125 


7400 


0.141 


0.116 


0.022 


0.024 


0.019 


0.017 


0.093 


0.073 


0.070 


0.095 


0.105 


0.116 


7600 


0.105 


0.106 


0.016 


0.019 


0.014 


0.012 


0.058 


0.068 


0.062 


0.086 


0.095 


0.106 


7800 


0.077 


0.094 


0.011 


0.012 


0.008 


0.006 


0.038 


0.063 


0.054 


0.077 


0.084 


0.094 


8000 


0.070 


0.077 


0.009 


0.009 


0.006 


0.004 


0.015 


0.052 


0.052 


0.071 


0.077 


0.086 


8200 


0.060 


0.058 


0.011 


0.011 


0.009 


0.008 


-0.006 


0.043 


0.052 


0.066 


0.071 


0.080 


8400 


0.050 


0.042 


0.011 


0.011 


0.011 


0.011 


-0.025 


0.051 


0.051 


0.063 


0.066 


0.074 


8600 


0.035 


0.032 


0.008 


0.008 


0.008 


0.009 


-0.035 


0.043 


0.042 


0.051 


0.053 


0.060 


8800 


0.017 


0.025 


0.002 


0.002 


0.003 


0.003 


-0.042 


0.024 


0.027 


0.035 


0.035 


0.040 


9000 


0.003 


0.019 


-0.002 


-0.003 


-0.002 


-0.002 


-0.060 


0.026 


0.021 


0.027 


0.026 


0.029 


9200 


0.002 


0.010 


-0.003 


-0.004 


-0.003 


-0.003 


-0.048 


0.009 


0.014 


0.019 


0.016 


0.018 


9400 


-0.001 


0.004 


-0.002 


-0.003 


-0.003 


-0.003 


-0.051 


0.003 


0.012 


0.016 


0.012 


0.013 


9600 


-0.007 


-0.001 


-0.001 


-0.003 


-0.003 


-0.003 


-0.060 


0.002 


0.012 


0.014 


0.010 


0.009 


9800 


-0.016 


-0.005 


-0.001 


-0.003 


-0.004 


-0.004 


-0.069 


0.002 


0.012 


0.013 


0.008 


0.007 


10000 


-0.026 


-0.010 


-0.001 


-0.003 


-0.004 


-0.005 


-0.080 


0.000 


0.013 


0.014 


0.008 


0.006 



al. [1998 , 2002) and limited the parameter space to T e fj = 
2 000-10 000 K, logg = 0, 3 and 4.5 and solar metal abun- 
dance. We use a AO V model with T eff = 9500 K, logg = 4.0 
and [Fe/H] = 0.0 to define zero-colours between all filters. 

In Figures [6] and [7] we have plotted all filter curves used 
for the colour calculations, together with synthetic stellar at- 
mospheres of different temperatures (log g = 4.5). In Table [4] 
we summarize these calculations. 

3.4.2. Finding excess stars from Spitzer and ISO data 

The interstellar extinction is certainly less in the IRAC bands 
than in the near-IR, making it possible to separate normal and 
excess stars without first calculating the extinction to each 
source. In Fig.|8]we have used a Spitzer-only (4.5 and 8.0^m) 



colour-magnitude diagram to find YSO candidates, taking into 
account the exclusion of both non-excess stellar and extragalac- 
tic sources. This method has been shown to be very efficient 
(see e.g. J0rgensen et al. 2006 and Harvey et al. 2006). 

The first criterion assumes that the colour index [4.5] - [8.0] 
is insensitive to interstellar extinction. Thus, stars with [4.5] - 
[8.0] less than the largest colour index for normal stars (we get 
0.27 for giants and 0.20 for MS stars) plus some margin, we 
chose 0.5, are assumed to have no intrinsic IR excess. This also 
means that we identify all the stars with [4.5] - [8.0] > 0.5 as 
stars with intrinsic IR excesses and therefore YSO candidates. 

The second criterion, [8.0] > 14 - ([4.5] - [8.0]), suggests 
possible remaining extragalactic sources in the sample (sources 
that are obvious galaxies in our / or K$ mosaics have already 
been excluded). As can be seen in Fig. [8] five of our sources 
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Fig. 8. Colour-magnitude diagram using only Spitzer photometry (4.5 and 8.0^m). Sources marked with an asterisk show clear 
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background source in our Ks image and is located in the part of the diagram expected for a galaxy. The arrow indicates the 
expected interstellar extinction vector (Ak = 5 corresponds to Ay = 47.4). YSO candidates without excess are marked with a red 
square. 



are located below this (dashed) line. Interstellar extinction, al- 
though low at 8.0 j-im, may have pushed these sources below the 
line. Nevertheless, they must be considered less reliable YSO 
candidates than the others. 

Due to the sensitivity limit of the 8.0 /im observations 
(high contrast background from PAH emission bands that ef- 
fectively masks faint point sources at these wavelengths), there 
are sources in our list that lack detection at 8.0 /mi, but still may 
exhibit intrinsic mid-IR excesses. 

We find an excess at 8.0/mi for exactly 2/3 of the sources 
detected at 4.5 and 8.0 /im. It is however also noted that several 
of the non-excess sources are in fact YSOs, as confirmed by 
our spectroscopy. 



ner of Fig 12 stands out as being the brightest and most red 
source. This is the newly discovered outflow source (Galfalk & 
Olofsson 2007) that could not be shown to have excess in Fig. [8] 
because of saturation in the Spitzer data. 



3.4.3. Defining the interstellar reddening and 
identifying additional YSOs 

We first consider the I-J\sJ-Ks diagram (Fig . |9]l . Red aster- 
isks mark sources with a clear excess at 8 /mi, and green circles 
denote sources not detected in the [8.0] band. In this Figure, 
both the / and J band sensitivity limits the sample we can use. 



Figure 12 shows our ISOCAM colour-magnitude diagram This is shown in Figure 10 where the observed / and / band 



for all sources detected at 6.7 and 14.3 /mi. The symbols are 
the same as in the Spitzer plots. There is a very clear separation 
between non-excess and excess sources at these longer wave- 
lengths and an excellent agreement between excess sources at 
8.0 and 14.3//m. Source L1641N-172, in the upper right cor- 



fluxes have been calculated for increasing distances towards the 
cloud edge and increasing extinctions within the cloud. For late 
MS M-type stars both / and J limits the sample to a few magni- 
tudes of extinction in the J band. However, for earlier spectral 
types and younger PMS stars it is the / band that sets the limit 
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Fig. 9. Colour-colour diagrams using ground based and Spitzer photometry. The solid black curve represents model stellar atmo- 
spheres (logg = 4.5, Tgff = 2 500 - 10 000 K) and the brown lines the reddening band due to interstellar extinction. The slopes 
of the reddening vectors were found by fitting the 'blue' sources. Symbols are the same as in Fig. [8]with the addition of green 
circles that denote sources not detected at 8.0 fj.m and thus having unknown red/blue status in the [4.5]-[8.0] index. Error bars are 
shown for selected sources located on the red side of the reddening band (all such green circles with error bars are added to our 
list of YSO candidates). It is noted that very few of the sources show clear intrinsic excess at near-IR wavelengths (K s ). (Aj -5 
corresponds to A v = 19.4) 



12 M. Galfalk and G. Olofsson: A detailed study of the L1641N star formation region 



20 



\\ V V 


i i 1 i i i 1 i 


\ \ \ ^\ ^^^^ =:= 5a^ — -— ___K5 (0.1) 


_ 


\ \\\ ^^^^ a ===== = -_ K5 < 1 > 

- \ \ \ \^ ^~~~~ ;::: ^^^== === 48(MS) 


\ — 


" \ \v M5(0.1) 


\ \ 






K5(MS) 


\ \ _ 


MO (MS) 
^^^^^^^ 


\ \ \ — 
V\\ \v \ 


M5 (MS) 


\ \ \ \\. \ — 




D! 




T3 


-a 




O 

O 




I band limit 




■ i i I i i i I 





200 400 1 2 3 4 5 6 7 8 9 

Distance [pc] Absolute extinction A, [mag] 



6 



10 



14 



16 



18 





'X* i 


1 1 


I I 






~~~--~~^K5 (0.1) 








""~~~--^^M5(0.1) 
















V^~-^8(MS) 








M5(1) 








■~~~~~~~^K5 (MS) 








~~~~ — -~-~-__JVI0 (MS) 








M5(MS) 






2MASS limit 


d edge 




1 


1 


- 

Clou 





200 400 1 2 3 4 5 6 7 8 9 

Distance [pc] Absolute extinction A, [mag] 



Fig. 10. The completeness of the survey as determined by source absolute magnitude, distance and extinction in the / and / band 
observations. The age of the stars are shown in parenthesis following the spectral type. The horizontal axis represents distance 
and J band extinction. Limiting magnitudes are marked by red dashed lines. (Aj = 5 corresponds to Ay = 19.4) 



M. Galfalk and G. Olofsson: A detailed study of the L1641N star formation region 



13 



h 1 1 1 1 1 r 



l 1 1 1 r 



~> — r 



-i 1 1 1 1 1 V 



3 — 
v> — 




145 
-* — i 



140 
— *— 



_l l l I l l l_ 



_l l l_ 



Ks - [4.5] 



sr-i r 



i r - 



-i — 7 r- 



— 




A . = 15 mag 



Ks - [8.0] 

Fig. 11. Additional colour-colour diagrams using ground based and Spitzer photometry. The symbols and curve/line styles are the 
same as in Figures [8] and [9] There is a clear separation between intrinsic non-excess and excess sources at both 4.5 and 8.0/v.m, 
with the non-excess sources agreeing well with the stellar model atmospheres and interstellar reddening. The blue sources are 
somewhat bluer that expected in the K s -[4.5] and K s -[3.6] indices (Fig. [9]), this is probably due to a small systematic shift in the 
Spitzer 3.6 and 4.5 fim magnitude calibrations. (Ay = 5 corresponds to Ay = 19.4) 



14 



M. Galfalk and G. Olofsson: A detailed study of the L1641N star formation region 



10 



32 
+ 



98 
+ 
+ 

168 



O 

172 



12 3 

[6.7] - [14.3] 

Fig. 12. Colour-magnitude diagram showing that most of the 
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the entire LI 64 IN. It is marked with a circle since it could not 
be shown to have excess in the Spitzer data because of satura- 
tion. 



to which stars can be used (because of the higher extinction at 
790 nm). 

The 'normal' stars in the / — / vs J — Ks diagram (blue 
plus signs) occupy a band defined by the reddening vector and 
the curve of model stellar atmospheres. The lower border of 
the reddening band is relatively well defined. The slope is 1.81, 
meaning that E(I- J) = 1.81E(J - K s ). It is interesting to note 
that only a few of the stars with excess at 8 fim (red asterisks in 
the Figure) have significant excess emission at Ks . 

We next turn to the J - Ks vs Ks - [3.6] diagram (Fig. [9] 
lower panel). Again, we can relatively well define the redden- 
ing vector from the normal stars, and we get E(J - Ks) - 
A.5E{K S - [3.6]). The reason for the 'blue' sources (plus signs) 
being slightly offset to the left of the reddening band (by about 
0.2-0.3 magnitudes) is probably due to the calibration of the 
3.6 /mi observations. We note that many of the stars with ex- 
cess emission at 8.0/im also show excess emission at 3.6 /jm. 
It is also clear that some of the objects not detected at 8 /jm 
have significant intrinsic excess emission at 3.6 fim. Taking the 
photometric accuracy for the individual sources into account, 
we can add these sources to the list of YSO candidates having 
intrinsic IR excess emission. The J - Ks vs Kg - [4.5] diagram 
(Fig. [TT) essentially confirms what is shown in the previous 
diagram. The slope of the reddening vector is 3.64 and thus 
E(J-Ks) = 3.64E(K s -[4.5]). 

Finally, we consider the Ks - [3.6] vs Ks - [8.0] diagram 



(Fig. 11 lower panel). Also in this case there is a well-defined 
reddening vector, E(K S - [3.6]) = 0.70E(K S - [8.0]). It is 
interesting to note that the stars with mid-IR excess emission 
are spread along a locus roughly parallel to the reddening vec- 
tor, indicating that these excess stars could well have intrinsic 
colours that differ far less than the observed colours. 



Sources not detected at 8.0 /mi that are shown to have in- 
trinsic mid-IR excess using colour-colour diagrams like these 
have been added to our list of YSO candidates. All sources with 
intrinsic mid-IR excess are marked by bold style source num- 
bers in the Tables. 



In Fig. 13 we summarize the reddening information derived 
above and as a comparison we show the theoretical extinction 
curve calculated by Draine et al. ( 120031 for R = A V /E(B-V) = 
5.5. The reason why the calculated extinction in the [3.6] and 
[4.5] bands are significantly lower than the observed ones is not 
quite clear, but at least part of the explanation must be absorp- 
tion bands of ices which are not included in the model. There 
may also be an additional dust component, not included in the 
model, which adds significantly to the mid-IR extinction, as a 
similar lack of a deep extinction minimum around 5-6 microns 
was observed by ISO in the diffuse interstellar medium towards 
the Galactic Centre (Lutz et al. U996b . 

We note that the ratio E(I - J)/E(J - K s ) = 1.81 corre- 
sponds to an exponential extinction law: 



1.25 



where 



1.81 = 



(A,-Aj) Aj/Aj-l (1.25/0.80/ -1 



(Aj-A Ks ) 1 -A Ks /Aj 1 - (1.25/2.14)0 
solving this equation gives /? = 1.58 and since 

E(I-J)=Aj-Aj=A^ 

we get an empirical absolute /-band extinction of 
Aj = 0.965£(/ - J) 



1 1.25 V* 




7 1.25 V* l 




-Aj=Aj 




\0.8oj 




AO.80] , 



(1) 



(2) 



(3) 



(4) 



Equation|4]is used to correct our / band photometry for ex- 
tinction, where E(I-J)q has been found from relating observed 
and model stellar colours in the I - J vs. J - Ks diagram. 

An extrapolation of this power law (J3 = 1.58) to zero wave 
number gives Ag s = 0.8E(J-Ks) (see Fig. 13 1. The theoretical 
value is similar, namely Ak s = 0.67 E(J - Ks). Even though 
there is a discrepancy between the observed and the theoretical 
extinction in the [3.6] and [4.5] bands, we adopt this latter value 



and we get the following extinction law 


A(I) 


= 3.48£(/-^) 


A(J) 


= 1.67 E(J-Ks) 


A(Ks) 


= 0.67E(J-K S ) 


A([3.6]) 


= 0.45£(/-^) 


A([4.5]) 


= OAOE(J-Ks) 


A([8.0]) 


= 0.35E(J-K S ) 



3.5. Photometric classification and corrections for 
extinction 

We return to the / - J vs J - Ks diagram (Fig. |9]l and note 
that, assuming no excess emission in these bands, it is possible 
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0.0 0.5 1.0 1.5 

MX [urn" 1 ] 

Fig. 13. Interstellar extinction curve from the optical (0.67 fim) through IR. The red plus signs represents measurements from our 
observations (/, J, K$, [3.6], [4.5] and [8.0]). The dotted line is a power law fit to the observed E(7 - J)/E(J - K$) and the solid 
line is a theoretical R = A V /E{B - V) = 5.5 extinction curve from Draine et al. (12003I I. 



to de-redden M-type stars and both determine their tempera- 
ture and the (absolute) extinction. As the spectral energy dis- 
tribution peaks in this wavelength region for these stars we can 
determine their luminosity with a reasonable accuracy. On the 
other hand, for hotter stars we need independent spectroscopy 
for the extinction correction. 

After de-reddening the M stars in the I — J vs / — Ks di- 
agram we obtain the intrinsic (7 - 7)o colours from where the 
de-reddening lines crosses the model curve. Using polynomial 
fits of bolometric correction BC(7) and T e ^- to intrinsic colours 
(/ - J)o, as shown in Figure 14 we then calculate the luminosi- 
ties and effective temperatures from photometry. The results 
are presented in Table [5] 

3.6. Spectral classification 

In total we have taken optical spectra (5780-8340 A) of 53 
sources in LI 64 IN, the observation log and results are pre- 
sented in Table [6] and the spectra themselves in Fig. 24 These 
stars were selected from our I band image and time constraints 
of our two spectroscopic observing runs, starting with spec- 
tra of the brightest I band sources and observing progressively 
fainter stars until the total exposure time needed for a single 
spectrum reached about 2 hours. 

In order to compare our observations to evolutionary mod- 
els, we need the effective temperature, T e ff- For this determi- 
nation we use synthetic spectra provided by P. Hauschildt (per- 
sonal communication) and in Table|6]we give the results, where 
we also include the photometrically determined temperatures. 



Table 5. Photometric mass and age estimates. 



The agreement is relatively good (see Fig 15 i and the scatter is 
about 200 K (1 cr). These GAIA-V2.0 LTE model spectra cover 
J eff = 2000 - 4000 K in steps of 100 K (late K and M spectral 
types). Smoothing was used to adapt the grid of high-resolution 
model spectra (0. 1 A resolution) to the resolution of the obser- 



No. 



A, 
(mag) 



*eff 

(K) 



m bo , 
(mag) 



L 

Oo) 



M 

CM©) 



Age 
(Myr) 



23 


1.69 


3560 


6.47 


0.21 


0.43 


3.14 


28 


1.37 


3110 


9.54 


0.012 


0.13 


18.71 


31 


2.79 


3550 


6.76 


0.16 


0.43 


4.98 


47* 


1.46 


3470 


8.78 


0.024 


0.40 


80.41 


64 b 


0.71 


3180 


9.96 


0.008 


0.16 


52.89 


80 


3.41 


3150 


6.74 


0.16 


0.19 


0.77 


84 


-0.41 












92 


2.11 


-2500 


8.68 


0.027 


< 0.1 


out 3 


118 


2.28 


3330 


5.76 


0.39 


0.25 


0.37 


128 


3.45 


3410 


5.36 


0.57 


0.28 


0.26 


133 


4.66 


3840 


4.03 


1.93 


0.43 


0.13 


138 a 


-4.61 












138 b 


-3.46 












140 b 


1.14 


3320 


9.74 


0.010 


0.26 


115.45 


141 


-1.40 












148 


5.05 


4000 


4.24 


1.59 


0.55 


0.31 


156 


-0.11 












157 


-0.97 












167 b 


3.11 


3750 


7.69 


0.07 


0.58 


52.56 


168 


4.37 


3630 


2.94 


5.25 


0.31 


0.01 


170 


1.87 


-2500 


9.00 


0.020 


< 0.1 


out" 


174 


2.07 


2940 


8.80 


0.024 


< 0.1 


-2" 


186 


4.12 


3810 


6.67 


0.17 


0.69 


12.21 



a Outside model grid (mass lower than 0. 10 Mq) 
b Close to MS 



vations. We assumed a surface gravity of log(g) = 4.5 for all 
stars and used the extinctions, Ay, obtained from photometry 
to deredden the spectra. 
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Table 6. Optical spectroscopy - observation log and results. 



No. 




(mag) 


i eff 


m/, / 
(mag) 


L 


M 


A op 
(Wiyr; 


vv Ha 

(M 
l A ) 


W, 


Pnni m pn f 

V^UllllllGllL 


5 


300 


~0 


6500 


2.47 


8.10 


1.61 


9.62 


4.13 




Early (F5) 


15 


1200 


0.20 


3900 (3640) 


5.07 


0.74 


0.56 


0.93 


-4.80 


0.58 




22 


300 


0.05 


5800 


1.75 


15.64 


2.59 


2.11 


2.15 




Early (G4) 


30 a 


600 


0.33 


3200 (3270) 


4.92 


0.85 


0.19 


0.03 


-9.66 






32 


300 


0.00 


6200 


2.85 


5.70 


1.54 


10.04 


3.07 




Early (F8) 


46 


1200 


0.00 


2900 (2510) 


6.53 


0.19 


0.11 


0.08 


-27.43 






47 d 


6000 


















Too faint 


55 


600 


0.03 


3700 (3420) 


5.23 


0.64 


0.42 


0.61 


-6.21 


0.40 




56 


3000 


0.24 


3200 (3130) 


8.16 


0.04 


0.20 


6.49 


-66.10 






60 d 


3600 


1.23 


4600 


6.36 


0.22 


0.75 


51.64 








61 


900 


0.01 


4100 


5.10 


0.72 


0.75 


1.93 


-2.29 


0.32 




64 d 


3600 


















Too faint 


65 


6000 


1.38 


3000 


7.07 


0.12 


0.13 


0.47 


-22.42 






66 a 


3600 


1.55 


4050 


4.85 


0.91 


0.66 


1.03 


-21.56 






66 b 


3600 


2.08 


3400 (3650) 


5.69 


0.42 


0.28 


0.44 


-56.56 






68 


3000 


0.28 


3050 (3190) 


8.26 


0.04 


0.13 


2.78 


-11.25 






73 


1200 


0.97 


3400 (3530) 


5.41 


0.54 


0.27 


0.27 


-39.32 






75 a 


2700 


1.39 


3750 (3900) 


3.86 


2.25 


0.37 


0.07 


-5.01 


0.20 




75 b 


1500 


1.80 


3750 (3640) 


4.70 


1.04 


0.41 


0.29 


-7.25 






78 


1200 


0.07 


3500 (3370) 


6.26 


0.25 


0.37 


1.73 


-3.44 


0.53 




84 


6000 


















Too faint 


98 


300 


~0 


5700 


2.71 


6.51 


2.02 


3.90 


2.03 




Early (G6) 


114 


1200 


0.18 


3800 


5.50 


0.50 


0.53 


1.41 


-44.80 


0.49 




117 


3600 


0.03 


3500 (3300) 


6.81 


0.15 


0.39 


4.31 


-5.27 






118 


1200 


















Too faint 


119 


1200 


0.20 


2900 (2830) 


8.53 


0.03 


<0.1 


out" 


-13.95 






120 


3600 


0.46 


3300 (3940) 


8.90 


0.02 


0.26 


30.92 


-10.39 






121 


1200 


1.26 


3200 (3960) 


5.65 


0.43 


0.20 


0.18 


-123.0 






142 


7200 


3.84 


3800 


3.51 


3.12 


0.39 


0.05 


-70.69 






144 c a 


1200 


0.53 


4000 


5.90 


0.34 


0.79 


6.08 


-18.29 


0.66 




144 b 


1200 


0.47 


4000 


6.70 


0.16 


0.76 


25.53 


-40.46 






145 


6000 


3.28 


3300 (3550) 


3.83 


2.32 


-0.20 


out" 


-12.67 






152 


6000 


1.85 


3800 


5.49 


0.50 


0.52 


1.38 


-243.34 


1.30 




160 


600 


0.68 


3600 (3250) 


6.08 


0.29 


0.43 


1.88 


-8.14 


0.88 




165 


3000 


















RefLNeb. 


168 


5400 


















Too faint 


171 


1200 


1.52 


3400 (3410) 


4.80 


0.94 


0.25 


0.09 


-16.46 






175 


5400 


2.18 


2650 (2940) 


6.39 


0.22 


<0.1 


out 1 








177 


1800 


0.06 


2900 


8.44 


0.03 


<0.1 


~1 


-13.70 


0.86 




179 


1800 


0.13 


3000 


7.58 


0.07 


0.12 


0.86 


-17.62 






184 


1200 


0.73 


2650 


7.34 


0.09 


<0.1 


out" 


-74.22 






185 


3000 


0.09 


3300 (3330) 


8.33 


0.04 


0.26 


14.71 








15/ 


1ZUU 


U.jy 




A OA 


1 ft7 
1.6/ 


U.J J 


V.Zj 


T. 1 ft 
-J. 1 


n Af\ 

U.'+U 




188 


3000 


0.28 


3100 (3270) 


6.74 


0.16 


0.17 


0.60 


-16.91 


0.58 




190 d 


3000 


1.16 


3600 (3400) 


7.94 


0.05 


0.80 


42.52 








192 


1200 


0.58 


3100 (3220) 


5.18 


0.67 


0.17 


0.03 


-36.55 






195 


1800 


0.62 


3400 (3290) 


5.95 


0.33 


0.29 


0.69 


-21.94 


0.18 




204 


1200 


0.34 


4200 


5.25 


0.63 


0.90 


3.54 


-7.31 


0.50 




205 


1800 


0.24 


3000 (3050) 


8.65 


0.03 


0.10 


2.85 


-36.12 






208 


1800 


0.53 


3050 (3380) 


8.21 


0.04 


0.13 


2.57 


-36.96 






210 a 


1200 


0.24 


5000 


3.99 


2.00 


1.55 


3.85 


-1.7 b 


0.30 




210 b 


2400 


0.06 


2800 


8.75 


0.03 


<0.1 


out" 


-3.12 






215 


1200 


0.09 


3100 (3020) 


7.11 


0.11 


0.17 


0.97 


-9.74 


1.06 





a Outside model grid (mass lower than 0. 10 M0) 

b P Cygni profile 

c Double source not resolved in J or K s 

d Close to MS 
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(i-J)o 



Fig. 14. Polynomial fits (solid curves) of bolometric correction BC( J) and effective temperature T e fj to intrinsic colour index 
(I—J)q. Model atmospheres (BaSeL v2.2) with a surface gravity of log g = 4.5 have been used for calculations in the temperature 
range 2 000-10 000 K with a step size of 100 K (circles). The blue curves are fits made using all effective temperatures, while the 
red curves are more accurate fits in the 2500-4000 K region (which we use in our photometric temperature determinations). 



For the temperature fits, we used TiO (Valenti et al. 1998) 
absorption band heads (spectral types earlier than M6), VO 
absorption features (later than M6) and the general shape of 
the spectra. If a photometrically obtained extinction differed 
markedly from that observed in a spectrum, we re-fitted Ay 
spectroscopically using the same extinction law as in the pho- 
tometry. In the visual comparison, we normalized the flux 
between the models and the spectra using a plateau around 
7500 A. The fit in temperature (and thus the spectral type) 
has roughly the same accuracy as the model steps themselves 
(100 K). Our effective temperatures agree fairly well with those 
also found from K band spectra in Hodapp & Deane ( 1993), es- 
pecially for early spectral types. A few of the late-type stars in 
our sample (L1641N-73, 145 and 192) have much lower tem- 
peratures (about 1000 K) than those found from the K band 
spectra. 

The luminosity of the spectroscopic sample is calculated in 
the same way as for the photometric sample, by de-reddening 
the stars in the I — J vs J — Ks diagram (yielding the extinc- 
tion Ay and intrinsic colour (/— J)q). However, if the extinction 
found photometrically does not fit the spectrum, Ay is instead 
estimated from a spectroscopic fit. Then follows calculations 
of the bolometric J band correction, BC(7), and finally the lu- 
minosity is calculated (assuming a distance of ~450pc). 



Most of the stars clearly have late type (mostly M) spec- 
tra typical of young dwarfs. One way to confirm their youth, 
and thus that they belong to the cluster, is the presence of 
Ha in emission (tracing accretion) and Li 1/16707 in absorp- 
tion (lithium has been destroyed in more evolved pre-main se- 



quence stars). In Fig. 25 we show some cases that have both, 



and in Table [6] we give the equivalent widths of these lines 
where applicable. Even though the spectral resolution of our 
spectra (R=830) is marginal for detecting the Li 6707 line, in 
particular as the S/N is in general poor due to the faintness of 
the stars, it was found present in many stars. 



With the exception of the four brightest stars, almost all of 
the spectra show Ha in emission and are of M-type. Since also 
main sequence M dwarfs frequently show Ha in emission, this 
by itself is not enough to classify them as YSOs. On the other 
hand, the combination of intrinsic IR excess and Ha emission 
is a reliable indication of YSO status. In addition, if the star 
exhibits interstellar reddening, it must be at least at the distance 
of the cloud, and we can judge from its luminosity if it is a PMS 
star. 




1000 



500 



8. 



-500 



-1000 
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Fig. 15. Comparison of spectroscopic and photometrically derived effective temperatures. The solid line indicates where 
T e fj(spec) = T e ff(phot). As can be seen in the lower panel, the residuals are fairly well centred on the solid line with a stan- 
dard deviation of a few hundred Kelvins. The red asterisks represent two sources, L1641N-120 and 121, having larger residuals 
than the other stars. This is probably caused by a small amount of intrinsic excess emission in the K s band (enough to keep them 
within the reddening band in Fig.|9]but yielding higher photometric temperatures). 



3.7. Contamination and multiplicity in the YSO sample 
3.7.1. Background sources 

An optical (or / band) image of L1641N shows that the extinc- 
tion is high across the whole cluster. Our / band image shows 
a total of 87 sources down to / ~22.5 mag while the K s im- 
ages show most of the Spitzer sources and many additional 
sources (limiting magnitude K$ ~ 19.5 mag). It is very likely 
that a number of Ks sources in fact are extragalactic objects. 
However, as can be seen in the comment column in Table [10] 
none of the / band sources appear to be extragalactic. As we 
need / band photometry (or spectra at similar wavelengths) to 
calculate the luminosities and effective temperatures used in 
the H-R diagram, most likely we have no extragalactic objects 
in the sample we used to calculate the MR 

The situation is similar for the [4.5] and [8.0] observa- 
tions, which are the main filters we use to find YSO candi- 
dates through intrinsic IR excess. The extinction is much less 
at these wavelengths. Even so, given the very high background 



(and high contrast) in these filters due to PAH emission (Fig. [5]), 
the sensitivity limit is severely affected. It is therefore unlikely 
that any extragalactic sources are included as YSO candidates. 
This is supported by the [8.0] vs [4.5] -[8.0] diagram (Fig.[8]l in 
which just a few sources are seen below the 'extragalactic con- 
tamination' line. In addition, all these sources have a point-like 
appearance in the K s band. 



Background giants is another possible source of contami- 
nation. However, for all YSO candidates with intrinsic IR ex- 
cess (Table [TTj i and for most of the non-excess YSO candi- 
dates (Table 12 1 we have other evidence for youth and cluster 
membership. Even though we cannot rule out the possibility of 
background giants among the non-excess YSO candidates, for 
which we lack spectra and for which the extinction Ay is high, 
these should be very few. 
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Table 7. The six earliest stars in our spectroscopic sample. 



No. 


T eff 


MS Type 


M(V) a 


V-J b 


M(J) 


m(J) 


A, 


Distance modulus 


MS distance c 


Comment 








(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


M 


(pc) 




5 


6500 


F5 


3.68 


0.80 


2.88 


9.78 


0.00 


6.90 


240 


Foreground star 


22 


5800 


G4 


4.70 


1.15 


3.55 


8.89 


0.05 


5.29 


114 


Foreground star 


32 


6200 


F8 


4.00 


1.00 


3.00 


9.97 


0.00 


6.97 


248 


Foreground star 


60 


4600 


K4 


6.92 


1.95 


4.97 


14.30 


1.23 


8.10 


417 


Belongs to the cluster 


98 


5700 


G6 


4.90 


1.18 


3.72 


9.76 


0.00 


6.04 


161 


Foreground star 


210 a 


5000 


Kl 


6.13 


1.57 


4.56 


11.06 


0.24 


6.26 


179 d 


Belongs to the cluster 



a From Table 15.7 in Allen's astrophysical quantities (20001. 

b From Kenyon & Hartmann ( 1995). 

c Luminosity distance; assuming that all the stars have reached the MS. 

d This is a young pre-MS star and thus much brighter than assumed for the MS luminosity distance. 



3.7.2. Early type and foreground stars 

Given the distance to L1641N (d~450pc), and the limiting 
magnitudes of our I, J and K s observations, we expect that all 



MS foreground stars are seen in the observations (see Fig. 10 1 
. Cluster membership status for the late type stars in our sam- 
ple is found through a combination of intrinsic IR excess (from 
dust), spectral lines (Ha in emission, Li 6707 in absorption) 
and extinction. This method is less applicable for earlier, more 
massive young stars (which have a much quicker star formation 
process). 

The earliest stars in our sample (F and G stars) have vir- 
tually no extinction and are much brighter than the rest of 
the stars (predominantly M-type stars). Since there are no 
Hipparcos data for these sources, we could not make distance 
estimates from parallaxes. 

The six earliest (and therefore among the most luminous) 
stars in our spectroscopic sample are listed in Table [7] Their 
high apparent brightness could be caused by a combination of 
nearby distance, low extinction and high intrinsic luminosity. 
Using spectra to calculate their extinctions and distances, clus- 
ter membership status can be investigated. We assume in the 
calculations that these stars have all reached the MS (because 
of their relatively high temperatures and thus higher masses), 
however, the coolest of these are more likely to yield closer 
MS luminosity distances than their actual distances since they 
are probably still pre-MS stars and thus brighter than assumed. 

Both L1641N-60 and 210A are YSO candidates, with mea- 
sureable J band extinction, and thus belong to the cluster even 
though the MS luminosity distance of 210A is only 179 pc (a 
bright PMS star). We will return to this later, as it is supported 
by their positions in the H-R diagram. 

From the Wainscoat et al. ([T992) model of the galaxy, we 
should expect about four foreground stars for spectral types 
F0-K5. This agrees with the four F-G stars we see in the sam- 
ple. However, we should then also expect roughly nine M-type 
foreground stars, whereas we see virtually no late-type fore- 
ground stars. These would show up as M stars with almost no 
extinction, but for most of these stars we instead have evidence 
for them being YSOs. For the four early type stars, we find no 



signature of youth. From this result, combined with the lack of 
extinction, we must assume that these F and G stars are fore- 
ground stars. 

3.7.3. Double sources 

Using our /, J, Ks and L' images we have visually detected 
eight double sources in the sample. These are listed in Table [8] 
including coordinates, separations and component photometry. 
We do not claim that these are physically bound binaries, al- 
though the ones with small separations most probably are. The 
criterion we use is that they are not resolved in the Spitzer ob- 
servations. 

Most of the doubles with fairly small separations were 
found by subtracting a fitted PSF from all sources in the sam- 
ple, sometimes revealing a second component. In Figure 16 
some of these subtractions are shown for both double sources 
(primary subtracted) and reflection nebulae (illuminating 
source removed to decrease contrast). 

In all Tables we have used the suffix a (primary) and b (sec- 
ondary) when we denote the components of a double source, 
and no suffix when a double is treated as a single source. We 
find a double frequency of ~9 % using all our YSO candidates, 
however, since we can only detect fairly bright and wide dou- 
bles (separations <: 076), this means that a lot of the multi- 
ple YSO systems are most likely not resolved. Also, the very 
young, deeply embedded (Class 0) sources are not detected in 
the /, J or Ks bands, which we use for our (wide) double source 
statistics. 

There is clear evidence in literature that the binary fre- 
quency decreases with stellar mass. Duquennoy et al. d 1 99 1 1 
found a binary frequency of ~57 % for G-dwarf primaries and 
Marchal et al. d20031 > a binary frequency of -27 % for M0- 
M6 dwarfs. Ahmic et al. d20071 > report a binary frequency of 
~ll% for young very low mass objects (~0.03-0.15 M0) in 
the Chamaeleon I star formation region using adaptive optics 
on the VLT, confirming the trend for a lower binary frequency 
with decreasing mass. This agrees well with the study by Bouy 
et al. d20061 l that uses HST to look for BD binaries in the 
Pleiades (-13.3%). 
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Table 8. All known double 3 sources in L 164 IN. 



No. 


RA 


Dec 


I 




J 






Separation 


Comment 




(2000) 


(2000) 


(mag) 


(mag) 


(mag) 






30 a 


05:36:06.943 


-06:18:53.23 


13.383 ± 


0.009 


11.500 ± 


0.023 b 


10.496 ± 


0.002 


6.19 


High contrast pair 


D 


ID:3o:U6.j9s 


-U6:18:4y./o 


21.694 ± 


0.242 






16.637 ± 


0.170 






32 a 


05:36:06.870 


-06:23:34.94 


11.010 ± 


0.1 10 b 


09.965 ± 


0.021 b 


09.787 ± 


0.001 


1.93 




D 


U5:36:U6.991 


-1)0:23: 35.63 










1 1 1 C 1 i 

13. lot ± 


U.U30 






66 a 


05:36:11.474 


-06:22:21.82 


16.078 + 


0.010 


12.931 ± 


0.004 


10.761 ± 


0.010 


3.31 




D 


LD. jo. 1 1.440 


-U0.zZ.2D.lU 


17.386 + 


0.011 


14.121 ± 


0.020 


12.123 ± 


0.040 






75 a 


05:36:12.969 


-06:23:32.13 


14.139 ± 


0.025 


1 1.762 ± 


0.002 


10.181 ± 


0.010 


2.78 




b 


05:36:12.984 


-06:23:29.35 


15.980 ± 


0.060 


13.001 ± 


0.040 


11.162 ± 


0.010 






110 a 


05:36:17.859 


-06:22:28.55 










15.650 + 


0.070 


0.65 




b 


05:36:17.906 


-06:22:28.07 










16.399 ± 


0.100 






138 a 


05:36:21.044 


-06:21:53.88 


19.004 ± 


0.050 


13.749 ± 


0.010 


10.315 ± 


0.010 


1.61 


observed in L' band 


b 


05:36:20.961 


-06:21:52.84 


19.934 ± 


0.100 


14.221 ± 


0.010 


11.443 ± 


0.010 






144 a 


05:36:21.880 


-06:26:01.71 


14.430 ± 


0.030 


12.267 ± 


0.026 b 


10.262 + 


0.028 b 


1.03 




b 


05:36:21.814 


-06:26:02.08 


15.107 ± 


0.150 














210 a 


05:36:32.384 


-06:19:19.85 


12.440 ± 


0.100 


1 1.059 ± 


0.024 


09.860 ± 


0.027 


5.27 


High contrast pair 


b 


05:36:32.687 


-06:19:22.54 


17.170 ± 


0.150 


14.973 ± 


0.117 


14.140 ± 


0.030 







a Unresolved with IRAC. 

b Total flux of both components. 



We thus expect a binary frequency in our YSO sample of 
20-30 %, suggesting that there are unresolved binaries in our 
sample. We do not detect any visually unresolved doubles in 
our spectroscopy, all stars can be fit with a single model - this 
does however not mean that these are all single stars. If the 
effective temperatures of both stars are roughly similar or the 
luminosity contrast is high, a double source could easily go 
unnoticed. 

3.8. The resulting YSO sample and H-R diagram 

The full list of YSO candidates is divided into two parts, YSOs 
with intrinsic IR excess (Table [TT| and YSOs that lack IR ex- 
cess (Table \T2$. In total, we find 89 YSO candidates. 

In addition, we have divided the sources into five cate- 
gories, depending on the main method for classifying them as 
YSO candidates. These are given in the Tables for each source. 
The first three categories are used for sources where we have 
spectra. Category I is for YSO candidates where Ha is seen in 
emission and Li I A61QH in absorption. Category II for Ha in 
emission together with IR excess. Category III are YSOs with 
non-zero extinction, implying that they are located inside the 
cloud. In addition, they are located on a PMS track in the H-R 
diagram when this extinction has been corrected for. 

For stars that are too faint for optical spectroscopy, we use 
category IV for sources with non-zero extinction, a point-like 
appearance in the /, J and K s images and a H-R position on a 
PMS track. And finally, category V are YSO candidates with 
intrinsic IR excess as the only evidence of youth. 

There are two additional YSOs in the list, L1641N-124 and 
172. Both have intrinsic IR excess but are also the suggested 
sources for the two giant H2 flows in L1641N. There are prob- 



ably many more YSOs in L1641N than in our YSO lists, since 
we have to detect each candidate in either the [5.8] and [8.0] fil- 
ters (for evidence of IR excess) or the / and J band filters (for 
extinction estimates). These Spitzer filters have a high contrast 
background from PAH emission bands that effectively masks 
faint point sources at these wavelengths, and the / and J bands 
are limited by interstellar extinction even at fairly low extinc- 



tions (see Fig. 10 1. 



The resulting H-R diagram for all YSO candidates, for 
which extinctions could be calculated is presented in Figure 17 
Evolutionary tracks (blue curves) are shown for 0. 10-1.80 Mq 
in steps of 0.05 and for 1.80-3. 00 M© in steps of 0.10 (Dotter 
et al. 2007). Isochrones are marked by red dashed curves and 
the extinction (Ay) to which the survey is complete is illus- 
trated by orange curves. All YSO candidates are plotted with 
green asterisks, with the addition of parentheses for the four 
suggested (early type) foregrund stars and source numbers for 
double sources. 

Using this H-R diagram we have determined the mass and 
age of each YSO candidate (Tables [5] and [6]). The luminosity 
function, mass function, effective temperature and age distri- 



butions for this sample are presented in Section 3.10 



3.9. Uncertainties in the H-R diagram 

There are several complicated sources of uncertainty in the age 
and mass determinations using a H-R diagram. The two most 
important observational challenges are variability and multi- 
plicity. 

As mentioned previously, young stars are often variable be- 
cause of variable accretion onto the star, large star spots and 
variability in the extinction towards the star - colours could 
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110 (Sep. 0.65") 




Reflection nebulosity Binaries 



Fig. 16. K s band images of reflection nebulosity (left panel) and double sources (right panel) with and without the main source 
(illuminating source and primary, respectively) removed using a normalized PSF from an isolated star. Each image has a field 
size of 34'.'7 x 34'.'7. 



therefore be affected by this and observations should ideally 
be made in all filters within a very short time frame. For fairly 
bright sources we have used 2MASS photometry for J-K s 
colours since these are simultaneous observations. Our / band 
observations, however, were made at a different epoch, adding 
an uncertainty to the I-J colours. Also, our J and Ks band ob- 
servations were made at two separate epochs. 

It is generally believed that many young stars are binary 
systems (e.g. 79-86 % in sources driving giant HH flows, see 
Reipurth 2000), and by using all our observations at different 
wavelengths we have detected 9 % of the sources in our YSO 
sample to be double sources. If one could detect and include all 
true binaries in the H-R diagram this would of course increase 
the number of available sources, but also increase the stellar 
age estimates due to the components being fainter than the to- 
tal flux from a binary. Observationally, this would require very 
high resolution imaging using adaptive optics or interferome- 



try, since spectroscopy probably would fail because of spectral 
variability in these young sources. 

Since the components of a binary star are expected to have 
formed at roughly the same time (except for binaries formed in 
gravitational encounters between two pre-MS stars) we expect 
that the ages are similar for both components. As can be seen in 
the H-R diagram, although some of the double sources (espe- 
cially L1641N-66) have components on isochrones of similar 
age, there are double sources with components that do not ap- 
pear coeval. The fact that L1641N-210b is located at a lower 
mass than what is covered by the evolution tracks and the de- 
creased accuracy for older stars (such as L1641N-144) makes 
a comparison difficult. 

Another uncertainty in the H-R diagram is the distance 
module. The distance to the Orion GMC, and therefore Orion 
A which L1641N belong to, is poorly determined, in part due to 
its extension. E.g. Brown et al. ( 1994) finds a distance of 320 pc 
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Fig. 17. H-R diagram of the full (spectroscopic and photometric) sample for which we could obtain effective temperatures and 
extinctions. The evolution track models (blue curves with masses indicated) were obtained from the Dartmouth Stellar Evolution 
Database (Dotter et al. [2007 ). Isochrones are drawn with red dashed lines and their ages shown on the right side of the grid. 
Observations are marked using green asterisks, with parentheses for the suggested four early foreground stars and source numbers 
given for double sources. The orange curves and numbers illustrate the completeness of the survey across the grid in terms of the 
A j extinction. 



to the near edge of the Orion A & B clouds and 500 pc to the far 
edge. We have used the canonical distance of L1641N, which 
is 450 pc. A closer distance would mean that all sources are 
moved further down in the diagram. This would not affect the 
mass determinations very much since the evolution tracks are 
almost vertical in the H-R diagram, except for older or more 
massive stars. 

Regarding the models, neither increased luminosity due to 
accretion nor excess emission from possible disks are included. 
Reflected light in the atmospheres of the stars could also play 
a role. The inclination of an YSO with a disk could be another 
source of uncertainty in the luminosity estimate, since the ex- 
tinction correction assumes that the disk and interstellar extinc- 
tion laws are the same. 

Most of the uncertainties mentioned in this section affects 
the luminosity determination and therefore mainly the age es- 
timate of a star, while the mass estimate is much more accurate 



(especially for low mass stars) because of the general direction 
of the evolution tracks. The uncertainty in mass is particularly 
small for the spectroscopic sample where the effective temper- 
atures have been determined to within 100 K. 



3.10. The Mass Function 

We have used both our spectroscopic and photometric sam- 
ple to calculate the cluster mass function in L1641N. Stars for 
which we have spectra have more accurate estimates of stellar 
mass, given that T e fj is fit directly from the spectra. 



In Figure 18 we have plotted the resulting Luminosity 



Function (LF), effective temperature distribution, age distribu- 
tion and Mass Function (MF) for the full sample, and two MFs 
with different completeness in Ay, minimum mass and maxi- 
mum age (given in the Figure). From the T e fj distribution it is 
clear that most stars in the sample are M-type stars. 
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1 
log (Age) [Myr] 

Mass (Aj < 1, age < 25 Myr, M/Msun > 0.10) 



Mass (Aj < 2, age < 10 Myr, M/Msun >0.15) 




Fig. 18. Resulting temperature, luminosity, age and mass functions for the YSO sample shown in Fig. [171 The first four panels 
include all sources in the H-R diagram while the last two panels show the MF for the completeness limits given above the 
respective panels (the limiting mass is also illustrated by the red dotted lines). The red, hashed region in the age plot illustrates 
stars with intrinsic IR excess (disk stars). It is interesting to note that even though the excess YSOs are present at all ages, seven 
out of the eight youngest YSOs in the sample have IR excess. 



Looking at the age distribution, it is clear that intrinsic IR 
excess (disk) stars are present at all ages. We note that seven 
out of the eight youngest YSOs in the sample have IR excess. 
The median age of all YSOs is found to be ~ 1 Myr. 

The MF that includes all sources in the sample is biased 
by bright (young and/or hot) sources being seen at higher ex- 
tinctions than faint sources. The two additional MFs are un- 
biased within the completeness limits given in the last two pan- 
els in Fig. 



18 



and illustrated by the red dashed lines, mark- 
ing the minimum mass for which the MFs are complete. The 
mass distribution peaks at about 0.3-0.4 Mq for all three MFs, 
and is fairly flat for masses lower than the peak until the mass 



completeness limits are reached in both un-biased MFs. The 
mass function is thus broadly consistent with IMFs found else- 
where in literature, showing a turnover or plateau at a few times 
0.1M© (Elmegreen 120081 : Elmegreen et al. [2008l > - keeping 
in mind the limitations set by small number statistics, limited 
mass range, unresolved binaries, more low- than high-mass 
members leaving the cluster due to close encounters but also 
a possible counter-acting selection effect of low-mass YSOs 
being easier to find than higher mass YSOs due to varying disk 
lifetimes. 

The methods we use to find our YSO candidates are differ- 
ently sensitive to extinction (e.g. mid-IR excess detects YSOs 
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at much higher extinctions than optical spectroscopy) and the 
processes involved (e.g. YSOs that have dispersed their disks 
are harder to find due to no IR excess). But, we need /, J 
and Ks photometry for all sources that we include in the MF, 
since this is used to calculate the absolute extinction Ay and 
in turn the stellar luminosity (also the effective temperature 
for sources where we lack spectra). Therefore, both un-biased 
MFs are complete to the limits in A y, age and mass given in 



Fig. 18 These limits were found using the (orange) complete- 



ness curves in Fig. 17 which in turn where calculated using a 



grid of stellar temperatures and luminosities as input parame- 
ters (combined with the filter transmission curves and limiting 
magnitudes). The K s observations are deeper and the least sen- 
sitive to extinction of the three filters, therefore the / and J pho- 
tometry sets the completeness limits (as illustrated in Fig. |T0||. 

Given the large amount of sources confirmed as YSOs from 
our spectra, there is probably a large number of deeply embed- 
ded YSOs in the full L1641N source catalogue that we miss in 
our MF (since we need / and / photometry to include them). In 
other words, we cannot go much deeper than a few magnitudes 
in Ay if we want to include sufficiently low mass YSOs in the 
statistics. We also note that the number of sources in each mass 
bin is fairly small and that there may be a problem in general 
with cluster IMFs if the velocity dispersion of stars in a clus- 
ter (as will be discussed in the following section) depends on 
stellar mass. 



3.11. Star formation history 

From the H-R diagram (Fig.[TT| we can obtain the stellar ages 
and thus the star formation (SF) history of most stars in the H- 

Hin 



R sample. The SF history is shown at two scales in Figure 19 



a cumulative manner, the number of stars formed in the sample 
versus time. The time spans are 120 and 5 Myr with 100 bins, 
meaning bin sizes of 1 .2 and 0.05 Myr, respectively. The blue 
dashed curves represent simple theoretical models that assume 
a gravitationally unbound cluster with constant star formation 
rate (one star in 3.7 x 10 4 yr, based on the last 1 Myr of our 
sample) and stellar velocity dispersions of 1 and 3 km/s, re- 
spectively. For simplicity, all stars are assumed to form at the 
centre of the cluster in this model. 

Given that the uncertainty in the H-R age determination in- 
creases drastically with age, due to the slow luminosity change 
with time close to the main sequence, we cannot say with con- 
fidence that the first few stars in the plot really are older than 
~50 Myr. However, below this age the estimates are more cer- 
tain and we see an accelerating increase in the number of stars 
with time. It is interesting to note the large number of old stars 
(for Orion), we detect 11 stars older than 10 Myr. There are 
at least three possible reasons for this; the star formation pro- 
cess could be very slow to start, some older stars could have 
drifted into LI 64 IN or star formation has actively been going 
on for a very long time but the velocity dispersion has moved 
most of the older stars out of the cluster while the younger stars 
(<1 Myr) are still close to their birth sites. There is however an 
alternative, observational explanation for these old stars. Some 
of the YSOs with disks might be seen almost edge-on, so that 



mostly reflected light is detected, making them look fainter and 
thus older (as hinted by the fact that some of the oldest YSOs 
show IR excess). Evidence for old members in the Orion clus- 
ter, using both H-R diagrams and models of Li burning can be 
found in e.g. Palla et al. (120071 1 where several stars appear to 
have ages of 15-30 Myr. 



Figure 20 shows the spatial age distribution. Four deeply 
embedded (and presumably very young) sources, two of these 
being the outflow sources for the major outflows, are marked 
by asterisks. Different symbols and colours (red being the 
youngest and blue the oldest) are used to illustrate eight dif- 
ferent age groups. The arrows show the direction and out- 
flow sources of the two major bi-polar outflows in the region 
(Galfalk & Olofsson 2007). There is no very clear relation be- 
tween age and spatial distribution seen in this plot, but it is 
evident that the most embedded (and presumably youngest) 
sources, including the two outflow sources, are located close 
to the centre of the cluster. 



In Figure 21 we have divided the YSOs into two groups 



(younger and older than 1 .5 Myr, respectively) and plotted their 
relative cumulative radial distribution versus distance from the 
centre of L1641N (using the 2 mm continuum dust peak of 
Chen et al. (1995) as the centre point). The sources in the 
younger group are clearly found closer to central LI 64 IN than 
the older group. A spatial separation is also found if e.g. 0.5 
or 1.0 Myr is used to divide the two groups. The exact cen- 
tre of the cluster is difficult to pinpoint without a wide-field 
study (Fig. [20] suggests that the centre might be somewhat to 
the southeast of the IRAS source), thus the best we can do is 
to assume that the dust peak (very close to the IRAS source) is 
located at the exact centre. 

Madsen et al. (2002) have used Hipparcos data to calculate 
the velocity dispersion of stars in young clusters and associa- 
tions. They find typical velocity dispersions of roughly 1 km/s. 
For LI 64 IN, a 50 Myr star with such a velocity would (at most) 



have moved about 40 times the field of view in Fig. 20 and a 
1 Myr star at most about one field size (if the velocity is per- 
pendicular to the line of sight). 



As can be seen in Fig. 19 the star formation rate is roughly 
constant during the last 1 Myr (the youngest age bin includes 



the four very embedded sources marked by asterisks in Fig. 20 1. 



We estimate a star formation rate of one star in 3.7 x 10 4 yr. The 
high current star formation rate agrees very well with the pres- 
ence of two major outflows (the outflow phase has a duration 
of ~10 5 yr ; Reipurth et al. 2001 ) and the very high concen- 
tration of HH objects in L1641N. These results agree with the 
scenario suggested by Hodapp et al. ( 1993) of ongoing succes- 
sive formation of individual stars (although over a larger time 
span than the suggested 1 Myr) with no clear sign of a star burst 
event. 

Our results support the possibility that star formation in 
L1641N could have been going on at a fairly constant rate for a 
very long time, but only appears to have an increased rate with 
time and a current maximum because of velocity dispersion, 
making more and more YSOs leave the cluster as they age. The 



measured curve in Fig 19 has a shape similar to a simple model 
curve that assumes a constant star formation rate over a very 
long time with constant velocity dispersion (3 km/s) and an un- 
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Fig. 19. Cumulative star formation history plot for the sample where ages could be determined. The blue (dashed) curves are 
simple theoretical models assuming a constant star formation rate, an unbound cluster and velocity dispersions of 1 and 3 km/s, 
respectively. 



bound cluster. If there is a relation between stellar mass and the 
probability of a star to leak out of a cluster with time, e.g. low- 
mass stars being more easily thrown out due to close encoun- 
ters (see e.g. Delgado et al. 2003, Sterzik and Durisen 2003 ), 
this could affect cluster IMFs in general (unless the surveyed 
region is very large) in addition to the star formation history. 



4. Summary and conclusions 

In this detailed study of the L1641N cluster, a nearby star for- 
mation region in Orion, we have used space based (Spitzer and 
ISO 3.6-14.3 yum) and ground based (/, 7, K s and L') imaging, 
as well as optical spectroscopy (5780-8340 A) of a large num- 
ber of sources in an effort to find YSO candidates. YSO can- 
didate status is found through a combination of mid-IR excess 
(presence of a disk), Ha emission (accretion), Li non-depletion 
(sign of youth) and extinction (cluster membership). 

We have calculated colours of normal stars for the IRAC, 
/, J and Ks filters and used a / - J vs J - Ks diagram (and 
spectroscopy) to calculate extinctions, luminosities and effec- 
tive temperatures. We investigate statistical properties of this 
YSO sample through double sources, extinction, effective tem- 
peratures, age distribution, spatial distribution, star formation 
history and the Mass Function (MF) for different completeness 
limits. The following results were obtained: 

1 . We detect a total of 2 1 6 Spitzer and / band sources within 
the region of our ISOCAM survey. 

2. Using the observations and criteria mentioned above we 
have found 89 YSO candidates in the sample. 



3. Most of the spectra show M-type stars, with Ha strongly 
in emission, that belong to the cluster. Although the spectral 
resolution is somewhat low (R=830) we also detect Li 1/16707 
in absorption in many of the sources. 

4. The four brightest stars in the / band observations, which 
are also shown to have the earliest spectral types (F and G 
stars), are very likely foreground stars to L1641N. 

5. We find that the interstellar extinction is well fit by a 
power law (index = 1.58) in the optical and near-IR, although 
in the mid-IR (Spitzer) the observed extinctions are higher than 
expected from theory. 

6. Extended PAH emission is seen in high contrast all over 
L1641N, we suggest that the illuminating source is t Orionis 
(an 09 III star 30' to the north, distance -400 pc). 

7. Using Spitzer, Ks and deep 2.12yum H2 observations 
we suggest that IRAS 05338-0624 (the only IRAS source in 
L1641N) is in fact the combined flux of at least four sources 
(and possibly some reflection nebulosity). 

8. We find a double source frequency of ~9 % using all our 
YSO candidates, however, we can only detect fairly bright and 
wide doubles (separations ^ 076), meaning that a large number 
of double/multiple YSO systems are probably not resolved. 

9. The median age of our YSO sample is ~1 Myr. 

10. The unbiased MFs peak at about 0.3-0.4 M0 with an 
essentially flat distribution for lower masses. 

11. We find a star formation history with an accelerating 
increase in the number of stars with time, but note that this 
could be caused by migration of stars with time from the cluster 
and that a velocity dispersion of a few km/s (which is typical in 
star formation regions) is enough to account for this effect. We 
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Fig. 20. Spatial distribution of stellar ages. The two major bi-polar outflows are illustrated by arrows, with the corresponding 
outflow sources (and two additional deeply embedded sources, probably very young) marked by red asterisks. 



detect some very old stars (for Orion) of up to ~50 Myr (1 1 of 
the stars are older than 10 Myr). 

12. A more or less constant star formation rate is found dur- 
ing the last ~1 Myr and estimated to be one star in 3.7 X 10 4 yr. 
This agrees with the presence of two major bi-polar out- 
flows (the outflow phase has an expected duration of ~10 5 yr 
;Reipurthetal. l2TXTTT) . 

13. We find a spatial separation between older and younger 
YSOs (e.g. using a dividing age of 0.5, 1.0 or 1.5 Myr) from 
the distribution of the two groups. This shows that velocity dis- 
persion is an important factor to consider when investigating 
the star formation history and possibly also the IMF (hence we 
calculate the MF). In fact, a star with the typical velocity of 
~1 km/s could leave L1641N in a few Myr. 

14. These results agree with a star formation process that 
has been slow to start. Another possibility is a roughly constant 
star formation rate over a very long time scale, where the star 
formation rate only appears to increase with time and have a 
current maximum because of the velocity dispersion, making 
more and more YSOs leave the cluster as they age. If there is 
a relation between stellar mass and velocity disperson, cluster 
IMFs in general could be affected unless the surveyed region is 
very large. 
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Fig. 21. Cumulative radial distribution of YSOs. The relative number of sources versus distance from the 2 mm continuum dust 
peak (Chen et al. 1 19951 1 at the centre of L1641N. The younger group (0-1.5 Myr) is clearly more concentrated towards the centre 
while the older group (> 1.5 Myr) is more spread out (similar to a uniform distribution, dashed curve). This supports the age 
determinations and shows that the velocity dispersion of the cluster members is large enough that the older YSOs have had time 
to move away significantly from their birth sites. 
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Fig. 22. Deep and 2.12/im H2 colour composite of the central and south east region of L1641N. All visible sources in our 
sample are marked by a circle. The field shown has a size of 4'05 x 435. 
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Fig. 23. Deep and 2.12 fim H2 colour composite of the central and north west region of L1641N. All visible sources in our 
sample are marked by a circle. The field shown has a size of 4'40 x 4' 37. 
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Table 9. Photometry of all sources - space based 



1NO. 


KA 


Dec 


F3.6 


F4.5 


F5.8 


^6.7 


Fg.o 


F 14.3 




(2000) 


(2000) 


mjy 


mjy 


mjy 


mjy 


mjy 


mjy 


1 


05:36:02.90 


-06:26:10.7 


0.068 ± 0.023 


0.040 ±0.015 










2 


05:36:02.90 


-06:25:55.7 


0.110 ±0.012 


0.086 ±0.010 










3 


05:36:03.61 


-06:25:28.4 


0.073 ±0.016 


0.056 ± 0.008 










4 


05:36:04.19 


-06:24:55.1 


0.023 ± 0.004 


0.034 ± 0.003 










5 


05:36:04.22 


-06:20:44.0 


36.840 ± 0.076 


25.158 ±0.065 


16.222 ± 0.278 




8.036 ± 0.350 




6 


05:36:04.41 


-06:22:24.9 


0.062 ± 0.004 


0.061 ± 0.004 










7 
8 


05:36:04.56 
05:36:04.60 


-06:23:59.4 
-06:25:01.7 


0.070 ±0.010 
0.031 ±0.005 


0.053 ±0.010 
0.018 ±0.005 










9 


05:36:04.68 


-06:21:31.8 


0.068 ± 0.004 


0.042 ± 0.004 










10 


05:36:04.78 


-06:21:05.4 


0.082 ± 0.005 


0.057 ± 0.003 










11 


05:36:04.79 


-06:19:30.4 


0.127 ±0.020 


0.079 ± 0.010 










13 


05:36:05.01 


-06:24:01.7 


0.857 ± 0.057 


0.610 ± 0.024 


0.420 ±0.124 








14 


05:36:05.09 


-06:23:01.8 


0.062 ± 0.004 


0.035 ± 0.004 










15 


05:36:05.16 


-06:25:25.4 


14.263 ± 0.057 


9.762 ± 0.043 


6.504 ±0.169 


5.15 ±0.94 


3.137 ±0.303 




16 


05:36:05.44 


-06:18:15.1 


0.198 ±0.005 


0.151 ±0.003 


0.068 ± 0.021 








17 


05:36:05.53 


-06:23:42.9 


0.083 ± 0.004 


0.077 ±0.010 










18 


05:36:05.69 


-06:24:40.8 


0.069 ± 0.009 


0.075 ± 0.005 


0.054 ± 0.020 








19 


05:36:05.80 


-06:25:36.6 


0.070 ± 0.009 


0.037 ± 0.005 










20 


05:36:05.89 


-06:20:54.8 


0.135 ±0.014 


0.088 ± 0.006 










21 


05:36:05.99 


-06:20:25.9 


0.051 ±0.010 


0.057 ±0.010 










22 


05:36:06.05 


-06:19:38.7 


74.833 ± 0.239 


63.854 ± 0.214 


40.482 ±0.196 




21.628 ± 0.285 




23 


05:36:06.15 


-06:25:52.7 


3.979 ± 0.094 


4.045 ± 0.047 


3.770 ± 0.294 




4.016 ± 0.538 




24 


05:36:06.32 


-06:19:55.3 


0.119 ±0.033 


0.122 ±0.026 


0.143 ±0.039 




0.167 ±0.028 




25 


05:36:06.34 


-06:22:40.5 


0.070 ± 0.014 


0.058 ± 0.009 










26 


05:36:06.47 


-06:20:08.4 


0.057 ± 0.014 


0.055 ± 0.007 










27 


05:36:06.58 


-06:23:01.8 


0.037 ± 0.005 


0.014 ± 0.004 










28 


05:36:06.58 


-06:21:22.1 


0.174 ±0.028 


0.138 ±0.015 










29 


05:36:06.74 


-06:25:58.2 


0.141 ±0.049 


0.134 ±0.025 










30 


05:36:06.77 


-06:18:51.6 


21.982 ±0.055 


14.378 ± 0.045 


9.778 ± 0.289 


7.07 ± 1.54 


5.708 ± 0.842 




31 


05:36:06.83 


-06:24:22.8 


1.401 ±0.028 


0.967 ±0.010 


0.702 ± 0.095 




0.518 ± 0.326 




32 


05:36:06.87 


-06:23:35.2 


36.012 ±0.125 


23.849 ± 0.074 


14.810 ±0.174 


12.51 ± 1.14 


8.872 ± 0.363 


5.30: 


33 


05:36:06.88 


-06:25:42.7 


0.157 ±0.039 


0.108 ±0.022 










34 


05:36:07.33 


-06:24:59.2 


0.081 ±0.011 


0.064 ± 0.006 










35 


05:36:07.38 


-06:20:24.8 


0.090 ±0.013 


0.065 ± 0.008 










36 


05:36:07.55 


-06:19:52.3 


0.061 ±0.003 


0.030 ± 0.005 










37 


05:36:07.59 


-06:19:42.7 


0.052 ±0.012 


0.024 ± 0.007 










38 


05:36:07.60 


-06:24:40.5 


0.041 ±0.012 


0.058 ± 0.009 










39 


05:36:07.61 


-06:19:31.4 


0.050 ± 0.021 


0.040 ±0.018 










40 


05:36:07.74 


-06:22:44.5 


0.081 ±0.023 


0.053 ± 0.015 










41 


05:36:07.81 


-06:25:06.6 


0.028 ± 0.004 


0.024 ± 0.002 










42 


05:36:07.82 


-06:18:34.7 


0.017 ± 0.003 


0.015 ± 0.003 










43 


05:36:07.86 


-06:21:32.2 


0.067 ± 0.012 


0.068 ± 0.005 










44 


05:36:07.90 


-06:24:16.2 


0.120 ±0.005 


0.093 ± 0.005 










45 


05:36:07.91 


-06:22:19.3 


0.069 ± 0.008 


0.035 ± 0.004 










46 


05:36:08.33 


-06:24:37.9 


8.426 ± 0.036 


6.700 ± 0.029 


5.296 ± 0.056 


5.55 ± 1.63 


4.855 ± 0.091 


3.49 ± 2.35 


47 
48 


05:36:08.34 
05:36:08.70 


-06:18:25.9 
-06:19:13.4 


0.291 ±0.018 
0.070 ± 0.006 


0.181 ±0.017 
0.059 ± 0.004 


0.114 ±0.040 








49 


05:36:08.74 


-06:25:30.0 


0.036 ± 0.007 


0.032 ± 0.003 










50 


05:36:08.95 


-06:23:49.3 


0.049 ± 0.008 


0.048 ± 0.004 










51 


05:36:08.96 


-06:25:03.4 


0.060 ± 0.004 


0.049 ± 0.004 










52 


05:36:09.04 


-06:22:33.1 


0.145 ±0.003 


0.105 ±0.007 


0.039 ± 0.025 








53 


05:36:09.11 


-06:20:20.4 


0.123 ±0.015 


0.098 ± 0.008 










54 


05:36:09.21 


-06:26:24.6 


1.127 ±0.016 


0.881 ±0.011 


0.590 ± 0.053 




0.301 ±0.128 




55 


05:36:09.49 


-06:18:36.3 


13.814 ±0.024 


9.106 ±0.027 


6.318 ± 0.086 


4.10 ± 1.65 


3.178 ±0.163 




56 


05:36:09.52 


-06:24:33.7 


1.265 ±0.032 


1.128 ±0.020 


0.811 ± 0.051 




0.677 ± 0.075 




57 


05:36:09.70 


-06:25:11.5 


0.143 ±0.004 


0.131 ±0.003 


0.072 ± 0.022 








59 


05:36:09.86 


-06:20:14.8 


0.087 ± 0.005 


0.068 ± 0.006 










60 


05:36:10.19 


-06:18:54.8 


2.109 ±0.020 


1.357 ±0.018 


1.109 ±0.059 




1.058 ±0.168 




61 


05:36:10.44 


-06:20:01.5 


12.399 ± 0.062 


8.143 ±0.052 


5.755 ± 0.249 


6.27 ± 1.50 


2.895 ± 0.451 




62 


05:36:10.66 


-06:18:42.0 


0.117 ±0.015 


0.079 ± 0.014 










63 


05:36:11.02 


-06:20:39.8 


0.241 ±0.024 


0.221 ±0.013 










64 


05:36:11.09 


-06:18:48.7 


0.231 ±0.014 


0.151 ±0.014 


0.158 ± 0.075 








65 


05:36:11.11 


-06:18:14.1 


1.136 ±0.025 


0.795 ± 0.014 


0.348 ±0.1 17 




0.170 ±0.166 




66 


05:36:11.45 


-06:22:23.8 


42.688 ± 0.245 


42.300 ± 0.377 


36.000 ±0.971 


37.54 ± 1.71 


36.302 ± 1.589 


34.17 ±2.78 


67 


05:36:11.45 


-06:22:46.4 


0.044 ± 0.006 


0.029 ± 0.008 










68 


05:36:11.67 


-06:24:58.4 


1.230 ±0.014 


1.003 ±0.009 


0.792 ± 0.058 




1.077 ± 0.097 




69 


05:36:11.98 


-06:21:02.6 


0.050 ±0.010 


0.055 ± 0.004 










70 


05:36:12.23 


-06:19:17.7 


0.026 ± 0.005 


0.014 ± 0.003 










71 


05:36:12.48 


-06:19:52.3 


0.037 ±0.010 


0.031 ±0.008 










72 


05:36:12.60 


-06:18:05.7 


0.024 ± 0.005 


0.018 ±0.003 










73 


05:36:12.61 


-06:23:39.6 


15.379 ±0.019 


12.705 ±0.023 


10.984 ± 0.067 


15.38 ± 1.25 


10.443 ± 0.099 


12.81 ±2.00 


74 


05:36:12.95 


-06:21:33.8 


0.065 ± 0.014 


0.033 ± 0.006 










75 


05:36:12.97 


-06:23:31.8 


46.053 ± 0.045 


33.121 ±0.035 


27.002 ±0.091 


27.78 ± 0.82 


21.177 ±0.063 


22.07 ± 1.83 


76 


05:36:13.00 


-06:20:00.3 


0.130 ±0.007 


0.099 ± 0.005 


0.049 ± 0.013 








77 


05:36:13.04 


-06:22:56.0 


0.146 ±0.005 


0.132 ±0.007 


0.086 ± 0.039 








78 


05:36:13.15 


-06:25:41.1 


4.036 ± 0.097 


2.704 ±0.051 


1.815 ±0.256 


2.29 ± 1.17 


0.791 ± 0.291 




79 


05:36:13.17 


-06:19:29.6 


0.090 ± 0.007 


0.078 ± 0.012 










80 


05:36:13.42 


-06:20:45.6 


1.672 ±0.057 


1.202 ±0.034 


1.000 ±0.157 




0.461 ± 0.264 




81 


05:36:13.45 


-06:23:59.6 


0.058 ± 0.008 


0.054 ± 0.007 










82 


05:36:13.64 


-06:20:23.7 


0.065 ± 0.033 


0.055 ±0.015 
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Table 9. - continued from previous page 



No. 


RA 


Dec 


F* £ 
'3.6 


F , c 


^5.8 


Ft ^ 
r 6.7 








(2000) 


(2000) 


mjy 


mly 


mjy 


mly 


mly 


mly 


83 
84 


05:36:13.69 
05:36:13.70 


-06:20:39.3 
-06:18:02.7 


0.102 ±0.026 
0.242 ± 0.012 


0.080 ±0.015 
0.152 ±0.020 










85 
86 


05:36:13.72 
05:36:13.77 


-06:19:14.7 
-06:25:53.0 


0.083 ± 0.004 
0.020 ± 0.007 


0.067 ± 0.002 
0.016 ± 0.005 










87 
88 


05:36:13.90 
05:36:13.92 


-06:20:32.1 
-06:20:08.7 


0.186 ±0.026 
0.074 ± 0.006 


0.124 ±0.014 
0.082 ± 0.004 


0.086 ± 0.033 








89 


05:36:13.93 


-06:25:20.0 


0.119 ±0.003 


0.093 ± 0.003 


0.104 ±0.021 








90 


05:36:14.19 


-06:21:44.1 


0.055 ±0.013 


0.037 ± 0.007 










91 


05:36:14.24 


-06:19:50.7 


0.109 ±0.006 


0.087 ± 0.003 










92 


05:36:14.33 


-06:20:18.5 


0.707 ± 0.010 


0.495 ± 0.005 


0.259 ±0.1 17 








93 


05:36:14.47 


-06:21:05.6 


10.876 ±0.114 


8.426 ± 0.045 


5.868 ±0.139 


3.97 ± 1.74 


2.530 ± 0.466 




94 


05:36:14.49 


-06:22:58.7 


0.058 ± 0.006 


0.053 ± 0.009 










95 


05:36:14.73 


-06:18:44.2 


0.087 ±0.016 


0.055 ± 0.012 










96 


05:36:14.84 


-06:21:20.1 


0.031 ±0.005 


0.030 ± 0.007 










97 


05:36:14.94 


-06:18:17.8 


0.306 ±0.011 


0.243 ± 0.010 










98 


05:36:15.22 


-06:23:55.7 


42.809 ± 0.217 


28.743 ±0.138 


19.023 ± 0.239 


15.06 ±2.08 


10.530 ± 0.456 


2.48: 


99 


05:36:15.48 


-06:21:38.1 


0.036 ± 0.015 


0.032 ± 0.009 










100 


05:36:15.57 


-06:20:18.6 


0.170 ±0.024 


0.120 ±0.011 


0.069 ± 0.055 








101 


05:36:15.84 


-06:18:02.2 


0.049 ± 0.005 


0.031 ±0.003 










102 


05:36:15.88 


-06:21:36.3 


0.050 ±0.011 


0.033 ±0.010 










103 


05:36:15.94 


-06:18:13.9 


0.110 ±0.017 


0.072 ± 0.015 










104 


05:36:15.96 


-06:21:59.7 


0.184 ±0.010 


0.182 ±0.007 


0.137 ±0.052 




0.106 ±0.062 




105 


05:36:16.13 


-06:20:39.3 


0.171 ±0.010 


0.138 ±0.006 










106 


05:36:16.54 


-06:25:49.7 


0.029 ± 0.003 


0.024 ± 0.003 










107 


05:36:16.69 


-06:20:50.7 


0.052 ± 0.007 


0.034 ± 0.005 










108 


05:36:16.73 


-06:18:17.2 


0.051 ±0.006 


0.043 ± 0.005 










109 


05:36:17.64 


-06:22:49.5 


1.871 ±0.105 


2.237 ±0.125 


2.594 ±0.159 




2.624 ± 0.099 




110 


05:36:17.86 


-06:22:28.6 


4.657 ± 0.224 


11.478 ± 1.097 


18.250 ± 1.947 


47.98 ± 4.75 


19.020 ± 1.504 


33.26 ± 3.69 


111 


05:36:18.03 


-06:25:10.5 


0.161 ±0.010 


0.149 ±0.007 


0.116 ±0.030 








112 


05:36:18.24 


-06:19:49.6 


0.032 ± 0.007 


0.038 ± 0.009 










113 


05:36:18.43 


-06:25:54.3 


0.245 ±0.011 


0.371 ±0.012 


0.687 ± 0.040 




0.893 ± 0.063 




114 


05:36:18.48 


-06:20:38.7 


19.624 ±0.168 


16.121 ±0.106 


15.206 ± 0.524 


15.23 ± 1.29 


14.438 ± 1.602 


15.63 ±2.24 


115 


05:36:18.49 


-06:22:13.6 


2.939 ±0.181 


15.455 ± 1.604 


29.897 ± 2.453 




38.856 ± 1.986 




116 


05:36:18.81 


-06:22:10.6 


20.132 ±0.487 


67.687 ± 1.910 


98.332 ± 1.769 


239: 


82.674 ± 0.904 


230: 


117 


05:36:18.86 


-06:22:04.4 


5.832 ± 0.940 


7.204 ± 2.213 










118 


05:36:18.88 


-06:26:21.2 


5.197 ±0.027 


3.543 ± 0.030 


2.631 ± 0.107 




1.673 ±0.178 




119 


05:36:19.08 


-06:22:50.8 


1.191 ±0.321 


1.534 ±0.360 


0.854 ± 0.272 








121 


05:36:19.38 


-06:25:51.3 


11.388 ±0.040 


10.712 ±0.024 


10.232 ± 0.217 


10.92 ± 2.57 


10.711 ± 0.445 


9.43 ± 4.08 


122 


05:36:19.50 


-06:24:38.3 


3.830 ±0.010 


4.443 ± 0.026 


3.952 ± 0.079 


3.38 ± 0.77 


3.984 ±0.166 


4.53 ± 1.35 


123 


05:36:19.50 


-06:22:20.2 


10.874 ± 0.636 


38.388 ± 1.408 


56.925 ± 2.587 




62.449 ± 2.021 




124 


05:36:19.50 


-06:22:12.1 


23.768 ± 1.945 


99.861 ±5.316 


192.786 ± 5.995 


408: 


Sat. 


297: 


125 


05:36:19.60 


-06:19:28.7 


0.051 ±0.017 


0.035 ± 0.023 










126 


05:36:19.78 


-06:22:22.6 


4.398 ± 0.053 


7.660 ± 0.076 


11.423 ±0.590 




11.442 ±0.756 




127 


05:36:19.83 


-06:20:27.4 


0.169 ±0.018 


0.1 14 ±0.012 


0.039 ± 0.042 








128 


05:36:19.93 


-06:18:27.5 


5.107 ±0.062 


3.442 ± 0.029 


2.377 ±0.150 


1.15 ±0.96 


1.030 ±0.202 




129 


05:36:20.05 


-06:24:00.4 


0.233 ± 0.023 


0.190 ±0.026 


0.150 ±0.068 








130 


05:36:20.40 


-06:19:49.7 


0.046 ± 0.003 


0.021 ± 0.004 










131 


05:36:20.41 


-06:25:19.0 


0.346 ± 0.015 


0.421 ±0.016 


0.412 ±0.047 




0.145 ± 0.060 




132 


05:36:20.47 


-06:21:29.2 


0.219 ± 0.035 


0.270 ± 0.063 


0.468 ± 0.102 




0.524 ± 0.068 




133 


05:36:20.49 


-06:23:22.3 


14.260 ± 0.206 


15.268 ± 0.293 


13.925 ± 0.242 


13.18 ±2.21 


11.186 ±0.266 


13.93 ±5.58 


134 


05:36:20.58 


-06:22:16.5 


1.302 ±0.230 


1.365 ±0.334 


1.870 ±0.379 


5.96: 


2.257 ± 0.308 


4.07 ± 1.58 


135 


05:36:20.93 


-06:25:51.4 


0.112 ±0.014 


0.1 16 ±0.040 










136 


05:36:20.94 


-06:26:30.9 


0.094 ± 0.004 


0.063 ± 0.005 










137 


05:36:20.95 


-06:22:38.8 


0.696 ± 0.020 


0.906 ±0.031 


0.911 ±0.022 




0.668 ± 0.038 




138 


05:36:21.02 


-06:21:53.0 


57.037 ± 0.352 


50.140 ±0.274 


50.621 ± 0.783 


47.85 ± 2.70 


43.367 ± 0.994 


45.02 ±9.14 


139 


05:36:21.07 


-06:18:37.1 


0.069 ± 0.008 


0.055 ± 0.006 










140 


05:36:21.17 


-06:23:31.1 


1.413 ±0.435 


1.349 ±0.396 


1.214 ±0.253 




1.034 ±0.321 




141 


05:36:21.17 


-06:18:04.1 


0.215 ±0.017 


0.180 ±0.015 


0.209 ± 0.075 




0.199 ±0.100 




142 


05:36:21.57 


-06:22:52.6 


50.632 ± 0.456 


60.789 ± 0.628 


64.558 ± 0.822 


55.58 ± 5.08 


69.837 ± 0.847 


84.86 ± 8.48 


143 


05:36:21.77 


-06:19:07.1 


0.034 ± 0.005 


0.030 ± 0.005 










144 


05:36:21.85 


-06:26:01.9 


70.374 ±0.151 


72.869 ± 0.090 


73.239 ± 0.265 


61.45 ± 1.63 


72.330 ± 0.555 


46.32 ± 7.33 


145 


05:36:21.87 


-06:23:30.1 


65.521 ± 4.265 


95.882 ±7.811 


140.288 ± 1.637 


161.96 ± 6.41 


209.150 ± 1.847 


269.03 ± 11.11 


146 


05:36:21.93 


-06:20:31.1 


0.051 ±0.007 


0.043 ± 0.005 










147 


05:36:21.96 


-06:20:59.3 


1.251 ±0.041 


1.429 ±0.012 


1.142 ±0.156 




0.490 ±0.219 




148 


05:36:22.09 


-06:24:22.7 


7.139 ±0.156 


5.220 ±0.113 


3.621 ± 0.283 


3.02: 


1.602 ±0.265 




149 


05:36:22.33 


-06:19:19.1 


0.156 ±0.006 


0.115 ±0.006 


0.056 ±0.013 








150 


05:36:22.37 


-06:25:45.0 


0.049 ± 0.034 


0.037 ± 0.053 










151 


05:36:22.42 


-06:25:24.4 


0.167 ±0.033 


0.145 ±0.030 


0.074 ± 0.032 








152 


05:36:22.46 


-06:23:44.9 


15.212 ±0.231 


16.617 ± 0.297 


12.888 ± 0.494 


15.91 ±2.82 


12.672 ±0.707 


14.17 ±3.73 


153 


05:36:22.62 


-06:21:27.7 


6.331 ±0.310 


7.565 ± 0.392 


7.626 ± 0.361 


5.70 ± 2.67 


5.976 ±0.358 




154 


05:36:22.74 


-06:18:53.9 


0.103 ±0.004 


0.058 ± 0.006 










155 


05:36:22.82 


-06:18:18.0 


0.135 ±0.011 


0.089 ± 0.005 










156 


05:36:22.88 


-06:24:09.3 


0.147 ±0.014 


0.095 ± 0.021 


0.153 ± 0.058 








157 


05:36:22.97 


-06:20:23.8 


0.478 ± 0.068 


0.485 ± 0.032 


0.420 ± 0.255 








158 


05:36:23.12 


-06:21:13.1 


0.069 ± 0.004 


0.047 ± 0.027 










159 


05:36:23.16 


-06:19:07.8 


0.023 ± 0.005 


0.011 ±0.005 










160 


05:36:23.26 


-06:19:37.5 


3.519 ±0.080 


2.329 ± 0.046 


1.412 ±0.168 




0.535 ± 0.266 




161 


05:36:23.33 


-06:18:19.8 


0.052 ± 0.015 


0.026 ± 0.009 










162 


05:36:23.60 


-06:24:51.7 


54.875 ± 0.885 


99.752 ± 1.770 


168.575 ± 2.870 


159.25 ± 5.66 


177.656 ± 3.565 


182.41 ± 7.02 


163 


05:36:23.65 


-06:21:46.7 


0.061 ±0.017 


0.038 ±0.017 
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Table 9. - continued from previous page 



No. 


RA 


Dec 


F* £ 


F , c 


^5.8 


Ft -7 

" 6.7 








(2000) 


(2000) 


mjy 


mjy 


mjy 


mjy 


mjy 


mjy 


164 


05:36:23.67 


-06:20:42.5 


0.172 ±0.015 


0.167 ±0.010 










165 


05:36:23.77 


-06:23:11.8 


64.073 ± 0.348 


89.058 ±0.281 


104.205 ± 0.825 


89.58 ± 3.30 


1 10.468 ± 0.530 


120.53 ± 7.64 


166 


05:36:23.80 


-06:19:26.4 


0.147 ±0.022 


0.068 ± 0.005 


0.059 ± 0.040 








167 


05:36:23.81 


-06:18:54.0 


0.621 ±0.011 


0.462 ± 0.010 


0.373 ± 0.057 




0.278 ±0.114 




168 


05:36:24.00 


-06:25:27.2 


34.052 ± 0.837 


25.393 ± 0.829 


17.454 ±0.685 


13.30 ± 1.27 


9.096 ±0.931 


2.09: 


169 


05:36:24.08 


-06:24:55.4 


5.096 ± 0.353 


8.465 ± 0.482 


13.930 ± 1.457 


18.08: 


19.462 ±2.153 


43.80: 


170 


05:36:24.23 


-06:18:03.6 


0.302 ± 0.032 


0.257 ± 0.041 


0.097 ± 0.053 








171 


05:36:24.48 


-06:22:23.2 


16.636 ±0.185 


14.679 ± 0.042 


14.420 ± 0.832 


5.08: 


13.349 ± 1.269 




172 


05:36:24.61 


-06:22:41.7 


49.554 ±0.109 


Sat. 


Sat. 


319.35 ± 6.58 


Sat. 


964.32 ± 16.00 


173 


05:36:24.63 


-06:20:11.8 


0.139 ±0.023 


0.097 ± 0.015 










174 


05:36:24.72 


-06:21:54.2 


0.682 ±0.018 


0.739 ±0.019 


0.832 ±0.108 




0.544 ± 0.205 




175 


05:36:24.83 


-06:25:07.5 


3.902 ±0.125 


2.763 ± 0.234 


3.274 ± 0.405 


6.47: 


2.717 ±0.193 


5.92: 


176 


05:36:25.07 


-06:19:49.8 


0.078 ± 0.006 


0.082 ± 0.005 










177 


05:36:25.28 


-06:23:07.2 


1.006 ±0.065 


0.758 ±0.100 










178 


05:36:25.39 


-06:26:26.9 


0.119 ±0.009 


0.082 ± 0.008 










179 


05:36:25.41 


-06:24:31.2 


2.431 ±0.252 


1.695 ±0.278 


1.837 ±0.155 




1.431 ±0.454 




180 


05:36:25.55 


-06:19:16.2 


0.082 ± 0.005 


0.052 ± 0.006 










181 


05:36:25.85 


-06:23:31.4 


0.240 ±0.010 


0.245 ±0.011 


0.213 ± 0.022 




0.188 ±0.129 




182 


05:36:25.86 


-06:20:55.2 


0.066 ± 0.007 


0.052 ± 0.003 










183 


05:36:25.88 


-06:24:58.6 


69.525 ±0.513 


85.481 ±0.653 


122.902 ± 0.859 


116.14 ±4.10 


181.314 ±0.739 


266.63 ± 10.40 


184 


05:36:26.08 


-06:26:15.5 


2.528 ± 0.052 


1.979 ±0.040 


1.102 ±0.102 




0.487 ± 0.085 




185 


05:36:26.10 


-06:24:51.7 


2.414 ± 0.367 


2.589 ± 0.506 


1.484 ±0.436 








186 


05:36:26.16 


-06:20:14.8 


1.043 ±0.022 


0.751 ±0.022 


0.666 ± 0.041 








187 


05:36:26.71 


-06:26:29.0 


31.035 ±0.142 


19.722 ± 0.077 


13.223 ±0.129 


7.19 ±3.14 


7.021 ±0.194 




188 


05:36:26.83 


-06:24:57.5 


5.443 ± 0.202 


3.528 ±0.181 


2.668 ±0.186 


2.69 ± 1.36 


2.452 ± 0.333 




189 


05:36:26.93 


-06:18:17.5 


0.356 ± 0.007 


0.242 ± 0.005 


0.211 ±0.028 








190 


05:36:27.01 


-06:21:07.6 


0.407 ± 0.046 


0.262 ± 0.021 


0.241 ±0.189 








191 


05:36:27.36 


-06:25:53.3 


0.016 ± 0.004 


0.013 ± 0.002 










192 


05:36:27.72 


-06:23:12.3 


19.366 ± 0.033 


16.140 ±0.056 


14.792 ± 0.177 


17.17 ± 1.99 


16.464 ±0.851 


14.86 ± 3.86 


193 


05:36:27.72 


-06:26:03.5 


0.059 ± 0.008 


0.043 ± 0.007 










194 


05:36:27.75 


-06:21:57.3 


0.769 ± 0.048 


0.720 ± 0.029 


0.840 ± 0.222 




0.643 ± 0.326 




195 


05:36:27.89 


-06:25:36.0 


5.263 ±0.121 


3.740 ± 0.094 


2.561 ± 0.320 


3.06 ± 1.12 


2.345 ± 0.558 


5.03 ±1.14 


196 


05:36:27.96 


-06:26:44.6 


0.331 ±0.005 


0.252 ± 0.004 


0.178 ±0.012 








197 


05:36:28.35 


-06:22:30.6 


0.039 ± 0.008 


0.024 ± 0.009 










198 


05:36:28.43 


-06:23:04.8 


0.131 ±0.012 


0.101 ±0.013 










199 


05:36:28.66 


-06:24:08.7 


0.097 ±0.011 


0.091 ±0.009 










200 


05:36:29.09 


-06:19:12.7 


0.044 ± 0.007 


0.028 ± 0.003 










201 


05:36:29.71 


-06:22:52.9 


0.123 ±0.010 


0.083 ± 0.010 










202 


05:36:29.89 


-06:24:29.2 


0.156 ±0.033 


0.535 ± 0.098 


1.125 ± 0.202 


3.00 ± 1.45 


1.304 ±0.192 


3.37: 


203 


05:36:29.89 


-06:19:52.3 


0.019 ±0.010 


0.014 ± 0.004 










204 


05:36:30.10 


-06:23:10.2 


49.132 ±0.122 


46.072 ±0.131 


41.567 ±0.499 


35.70 ± 2.27 


44.223 ± 1.316 


65.43 ± 6.75 


205 


05:36:30.50 


-06:23:56.6 


0.882 ±0.051 


0.650 ± 0.025 


0.609 ± 0.203 




0.961 ± 0.626 




206 


05:36:31.14 


-06:20:24.1 


0.063 ±0.011 


0.060 ± 0.007 










207 


05:36:31.26 


-06:23:49.9 


0.034 ± 0.009 


0.032 ± 0.004 










208 


05:36:31.84 


-06:23:23.1 


1.133 ±0.048 


0.859 ± 0.038 


0.737 ± 0.073 




0.643 ± 0.131 




209 


05:36:32.14 


-06:19:55.4 


0.027 ± 0.006 


0.016 ± 0.004 










210 


05:36:32.52 


-06:19:21.2 


47.905 ±0.150 


38.587 ±0.117 


28.364 ±0.160 


31.62 ± 1.51 


25.201 ±0.184 


21.16 ±2.62 


211 


05:36:32.55 


-06:19:43.1 


0.043 ± 0.004 


0.038 ± 0.007 










212 


05:36:32.68 


-06:21:47.0 


0.148 ±0.005 


0.138 ±0.005 










213 


05:36:32.91 


-06:22:33.7 


0.130 ±0.003 


0.119 ±0.004 


0.082 ± 0.030 








214 


05:36:33.06 


-06:21:05.0 


0.187 ±0.009 


0.145 ±0.006 


0.127 ±0.061 








215 


05:36:33.81 


-06:19:27.0 


3.025 ± 0.039 


2.121 ±0.028 


1.381 ±0.044 




0.450 ± 0.075 




216 


05:36:33.91 


-06:20:57.8 


0.081 ±0.019 


0.064 ±0.011 
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Table 10. Photometry of all sources - ground based 



1NO. 


KA 


Dec 


/ 


./ 








(2000) 


(2000) 


mag 


mag 


mag 


Comment 


1 


05:36:02.90 


-06:26:10.7 


21.08 ±0.18 








5 


05:36:04.22 


-06:20:44.0 


Sat. 


9.78 ± 0.03 


9.60 ± 0.02 




7 
8 


05:36:04.56 
05:36:04.60 


-06:23:59.4 
-06:25:01.7 






17.07 ± 0.12 
17.55 ± 0.20 




12 


05:36:04.93 


-06:25:54.0 


22.64 ± 0.70 






Artefact? 


13 


05:36:05.01 


-06:24:01.7 


20.89 ±0.10 




14.24 ± 0.02 




14 


05:36:05.09 


-06:23:01.8 






17.11 ±0.27 




15 


05:36:05.16 


-06:25:25.4 


13.14 ±0.01 


11.82 ±0.03 


10.90 ± 0.01 




17 


05:36:05.53 


-06:23:42.9 






17.20 ±0.17 




18 


05:36:05.69 


-06:24:40.8 






17.62 ± 0.38 


Ext/Gal 


19 


05:36:05.80 


-06:25:36.6 






17.06 ± 0.16 




20 


05:36:05.89 


-06:20:54.8 


20.25 ± 0.06 




15.98 ± 0.04 




22 


05:36:06.05 


-06:19:38.7 


Sat. 


8.89 ± 0.03 


8.63 ± 0.02 




23 


05:36:06.15 


-06:25:52.7 


17.54 ± 0.03 


14.60 ± 0.04 


12.81 ±0.01 




24 


05:36:06.32 


-06:19:55.3 






16.24 ± 0.05 


Ext/Gal 


25 


05:36:06.34 


-06:22:40.5 






17.35 ±0.14 




26 


05:36:06.47 


-06:20:08.4 


22.98 ± 0.80 




16.78 ± 0.08 




27 


05:36:06.58 


-06:23:01.8 






17.40 ±0.20 




28 


05:36:06.58 


-06:21:22.1 


20.31 ±0.06 


17.13 ±0.22 


15.54 ± 0.02 




29 


05:36:06.74 


-06:25:58.2 


21.60 ±0.16 




16.35 ± 0.06 




30 


05:36:06.77 


-06:18:51.6 


13.38 ±0.01 


11.50 ±0.03 


10.49 ± 0.01 


Double 


31 


05:36:06.83 


-06:24:22.8 


20.07 ± 0.06 


15.98 ± 0.07 


13.56 ±0.01 




32 


05:36:06.87 


-06:23:35.2 


11.01 ±0.11 


9.97 ± 0.03 


9.76 ± 0.01 


Double 


33 


05:36:06.88 


-06:25:42.7 


21.44 ±0.26 




16.40 ± 0.05 




34 


05:36:07.33 


-06:24:59.2 


21.99 ± 0.32 




16.26 ± 0.07 




35 


05:36:07.38 


-06:20:24.8 


22.15 ±0.24 




16.55 ± 0.03 




36 


05:36:07.55 


-06:19:52.3 






18.15 ±0.14 




37 


05:36:07.59 


-06:19:42.7 






18.70 ±0.14 




38 


05:36:07.60 


-06:24:40.5 






18.43 ± 0.33 




39 


05:36:07.61 


-06:19:31.4 


22.92 ± 0.35 




16.72 ± 0.04 




40 


05:36:07.74 


-06:22:44.5 


22.19 ±0.38 




16.52 ±0.04 




43 


05:36:07.86 


-06:21:32.2 






16.74 ± 0.06 


Ext/Gal 


44 


05:36:07.90 


-06:24:16.2 






16.49 ± 0.08 




45 


05:36:07.91 


-06:22:19.3 






17.58 ± 0.07 


Ext/Gal 


46 


05:36:08.33 


-06:24:37.9 


15.07 ± 0.02 


12.70 ± 0.03 


11.72 ±0.01 




47 


05:36:08.34 


-06:18:25.9 


19.42 ± 0.06 


16.63 ±0.15 


14.97 ±0.11 




48 


05:36:08.70 


-06:19:13.4 






17.66 ± 0.05 




49 


05:36:08.74 


-06:25:30.0 






17.60 ±0.14 




50 


05:36:08.95 


-06:23:49.3 






17.61 ± 0.14 




51 


05:36:08.96 


-06:25:03.4 






17.43 ±0.12 




52 


05:36:09.04 


-06:22:33.1 






16. 10 ±0.02 




53 


05:36:09.11 


-06:20:20.4 


21.84 ±0.19 




15.93 ± 0.02 




54 


05:36:09.21 


-06:26:24.6 




16.73 ±0.13 


13.96 ±0.01 




55 


05:36:09.49 


-06:18:36.3 


13.11 ±0.03 


11.75 ±0.03 


10.91 ± 0.03 




56 


05:36:09.52 


-06:24:33.7 


16.62 ± 0.02 


14.64 ± 0.04 


13.70 ±0.01 




57 


05:36:09.70 


-06:25:11.5 






17.13 ±0.16 




58 


05:36:09.72 


-06:19:29.0 


22.47 ± 0.34 




17.97 ± 0.05 




59 


05:36:09.86 


-06:20:14.8 


21.06 ± 0.08 




16.22 ± 0.03 




60 


05:36:10.19 


-06:18:54.8 


16.48 ± 0.02 


14.30 ± 0.03 


12.90 ±0.01 




61 


05:36:10.44 


-06:20:01.5 


12.79 ±0.01 


11.66: 


10.85 ± 0.01 




62 


05:36:10.66 


-06:18:42.0 


22.59 ± 0.74 




17.01 ±0.12 




63 


05:36:11.02 


-06:20:39.8 






16.04 ± 0.02 


Ext/Gal 


64 


05:36:11.09 


-06:18:48.7 


19.32 ± 0.04 


16.91 ±0.17 


15.70 ± 0.20 




65 


05:36:11.11 


-06:18:14.1 


16.93 ± 0.02 


14.64 ± 0.04 


13.79 ± 0.04 




66 


05:36:11.45 


-06:22:23.8 


15.79 ± 0.02 


12.62 ± 0.01 


10.49 ± 0.02 


Double 


67 


05:36:11.45 


-06:22:46.4 






18.50 ±0.08 




68 


05:36:11.67 


-06:24:58.4 


16.69 ± 0.02 


14.74 ± 0.02 


13.77 ± 0.01 




69 


05:36:11.98 


-06:21:02.6 






17.05 ± 0.04 




70 


05:36:12.23 


-06:19:17.7 






17.79 ± 0.06 




71 


05:36:12.48 


-06:19:52.3 


21.46 ±0.11 




16.71 ±0.03 




73 


05:36:12.61 


-06:23:39.6 


14.96 ± 0.01 


12.74 ± 0.04 


11.36 ±0.03 




74 


05:36:12.95 


-06:21:33.8 






18.27 ± 0.14 




75 


05:36:12.97 


-06:23:31.8 


13.96 ±0.04 


11.46 ±0.02 


9.81 ± 0.01 


Double 


76 


05:36:13.00 


-06:20:00.3 






16.32 ± 0.03 


Ext/Gal 


77 


05:36:13.04 


-06:22:56.0 






16.80 ± 0.05 




78 


05:36:13.15 


-06:25:41.1 


14.58 ± 0.01 


13.11 ± 0.03 


12.25 ± 0.01 




79 


05:36:13.17 


-06:19:29.6 


20.80 ±0.10 




16.30 ± 0.02 




80 


05:36:13.42 


-06:20:45.6 


21.61 ±0.16 


16.37 ±0.12 


13.62 ±0.01 




81 


05:36:13.45 


-06:23:59.6 






17.79 ± 0.20 


Ext/Gal 


82 


05:36:13.64 


-06:20:23.7 






18.18 ±0.10 




83 


05:36:13.69 


-06:20:39.3 






16.78 ± 0.03 




84 


05:36:13.70 


-06:18:02.7 


19.20 ± 0.04 


16.57 ±0.13 


15.53 ±0.17 




85 


05:36:13.72 


-06:19:14.7 






17.93 ± 0.08 




86 


05:36:13.77 


-06:25:53.0 






17.98 ± 0.24 




87 
88 


05:36:13.90 
05:36:13.92 


-06:20:32.1 
-06:20:08.7 






15.83 ± 0.02 
18.19 ± 0.16 




89 


05:36:13.93 


-06:25:20.0 






17.00 ±0.11 




90 


05:36:14.19 


-06:21:44.1 






17.20 ±0.04 




91 


05:36:14.24 


-06:19:50.7 






17.01 ± 0.04 


Ext/Gal 


92 


05:36:14.33 


-06:20:18.5 


21.51 ±0.10 


16.93: 


14.75 ± 0.01 
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Table 10. - continued from previous page 



No. 


RA 


Dec 


j 


/ 








(2000) 


(2000) 


mag 


mag 


mag 


Comment 


93 


05:36:14.47 


-06:21:05.6 




15.79 ± 0.08 


11.77 ±0.01 




94 


05:36:14.49 


-06:22:58.7 






17.55 ±0.11 




95 


05:36:14.73 


-06:18:44.2 


21.99 ±0.29 




16.61 ± 0.06 




96 


05:36:14.84 


-06:21:20.1 






17.98 ± 0.07 




98 


05:36:15.22 


-06:23:55.7 


10.97 ±0.11 


9.76 ± 0.03 


9.42 ± 0.02 




99 


05:36:15.48 


-06:21:38.1 






19.19 ±0.28 




100 


05:36:15.57 


-06:20:18.6 






16.07 ± 0.02 




102 


05:36:15.88 


-06:21:36.3 






17.32 ± 0.08 




103 


05:36:15.94 


-06:18:13.9 


21.37 ±0.31 








104 


05:36:15.96 


-06:21:59.7 




17.55 ± 0.08 


16.12 ±0.04 




105 


05:36:16.13 


-06:20:39.3 






16.15 ±0.02 




106 


05:36:16.54 


-06:25:49.7 






17.84 ±0.11 




107 


05:36:16.69 


-06:20:50.7 






17.95 ± 0.08 




109 


05:36:17.64 


-06:22:49.5 






14.51 ±0.04 




110 


05:36:17.86 


-06:22:28.6 






15.28 ± 0.08 


Double 


111 


05:36:18.03 


-06:25:10.5 






16.50 ± 0.06 




112 


05:36:18.24 


-06:19:49.6 






18.08 ± 0.05 




113 


05:36:18.43 


-06:25:54.3 






17.15 ±0.12 




114 


05:36:18.48 


-06:20:38.7 


13.60 ± 0.01 


12.20 ± 0.03 


11.00 ±0.01 




117 


05:36:18.86 


-06:22:04.4 


15.16 ±0.02 


13.62 ± 0.01 


12.79 ± 0.02 




118 


05:36:18.88 


-06:26:21.2 


18.17 ±0.02 


14.35 ± 0.03 


12.23 ± 0.02 




119 


05:36:19.08 


-06:22:50.8 


17.14 ±0.02 


14.90 ± 0.03 


13.95 ±0.05 




120 


05:36:19.23 


-06:22:04.4 


17.58 ±0.03 


15.66 ± 0.05 


14.32 ±0.18 




121 


05:36:19.38 


-06:25:51.3 


15.37 ± 0.02 


13.16 ±0.03 


11.66 ±0.01 




122 


05:36:19.50 


-06:24:38.3 




17.49 ± 0.05 


13.45 ±0.01 




123 


05:36:19.50 


-06:22:20.2 






17.39 ±0.16 




124 


05:36:19.50 


-06:22:12.1 






16.40 ±0.15 




125 


05:36:19.60 


-06:19:28.7 






17.94 ± 0.07 


Ext/Gal 


126 


05:36:19.78 


-06:22:22.6 






15.33 ± 0.02 




127 


05:36:19.83 


-06:20:27.4 






16.22 ± 0.02 




128 


05:36:19.93 


-06:18:27.5 


20.09 ± 0.06 


15.17 ±0.04 


12.37 ± 0.03 




129 


05:36:20.05 


-06:24:00.4 






16.46 ± 0.03 




130 


05:36:20.40 


-06:19:49.7 






17.86 ± 0.05 




131 


05:36:20.41 


-06:25:19.0 






16.26 ± 0.09 




132 


05:36:20.47 


-06:21:29.2 






16.02 ± 0.03 




133 


05:36:20.49 


-06:23:22.3 


21.03 ±0.11 


15.22 ± 0.04 


11.77 ± 0.02 




134 


05:36:20.58 


-06:22:16.5 




17.22 ± 0.06 


14.23 ± 0.02 




135 


05:36:20.93 


-06:25:51.4 






16.68 ±0.10 




136 


05:36:20.94 


-06:26:30.9 






16.27 ± 0.09 




137 


05:36:20.95 


-06:22:38.8 




17.85 ± 0.09 


15.18 ±0.02 




138 


05:36:21.02 


-06:21:53.0 


18.62 ±0.16 


13.21 ±0.01 


9.99 ± 0.01 


Double 


139 


05:36:21.07 


-06:18:37.1 






17.02 ±0.15 




140 


05:36:21.17 


-06:23:31.1 


19.85 ± 0.04 


17.20 ± 0.07 


15.72 ± 0.04 




141 


05:36:21.17 


-06:18:04.1 


20.69 ±0.10 


16.98 ± 0.17 


15.38 ±0.15 




142 


05:36:21.57 


-06:22:52.6 


18.62 ±0.02 


13.87 ± 0.01 


10.88 ± 0.01 




143 


05:36:21.77 


-06:19:07.1 






18.00 ±0.08 




144 


05:36:21.85 


-06:26:01.9 


13.96 ±0.06 


12.27 ± 0.03 


10.26 ± 0.03 


Double 


145 


05:36:21.87 


-06:23:30.1 


17.99 ± 0.02 


13.40 ± 0.01 


10.71 ± 0.02 




146 


05:36:21.93 


-06:20:31.1 






17.12 ± 0.03 




147 


05:36:21.96 


-06:20:59.3 






15.35 ±0.01 




148 


05:36:22.09 


-06:24:22.7 


21.97 ±0.19 


15.88 ± 0.07 


12.16 ±0.03 




149 


05:36:22.33 


-06:19:19.1 






16.18 ±0.03 




150 


05:36:22.37 


-06:25:45.0 


22.42 ± 0.47 








152 


05:36:22.46 


-06:23:44.9 


17.15 ±0.02 


13.86 ±0.01 


11.70 ±0.02 




153 


05:36:22.62 


-06:21:27.7 






13.46 ± 0.02 




154 


05:36:22.74 


-06:18:53.9 






17.14 ±0.10 




156 


05:36:22.88 


-06:24:09.3 


21.01 ±0.08 


16.93 ± 0.05 


16.07 ± 0.02 




157 


05:36:22.97 


-06:20:23.8 


20.29 ± 0.06 


16.56 ±0.12 


15.21 ±0.01 




158 


05:36:23.12 


-06:21:13.1 






18.14 ±0.11 




159 


05:36:23.16 


-06:19:07.8 






18.59 ± 0.23 




160 


05:36:23.26 


-06:19:37.5 


14.83 ±0.01 


13.21 ±0.03 


12.36 ± 0.03 




162 


05:36:23.60 


-06:24:51.7 






13.50 ±0.06 




163 


05:36:23.65 


-06:21:46.7 




18.35 ±0.15 


17.23 ± 0.07 




164 


05:36:23.67 


-06:20:42.5 






16.68 ± 0.05 


Ext/Gal 


165 


05:36:23.77 


-06:23:11.8 


18.00 ± 0.06 


14.39 ± 0.07 


10.70 ±0.02 


Refl.Neb. 


166 


05:36:23.80 


-06:19:26.4 






16.31 ±0.05 




167 


05:36:23.81 


-06:18:54.0 


21.56 ±0.29 


17.31 ±0.22 


14.73 ± 0.02 




168 


05:36:24.00 


-06:25:27.2 


19.43 ± 0.06 


13.77 ± 0.04 


10.46 ± 0.03 




169 


05:36:24.08 


-06:24:55.4 






16.41 ± 0.09 




170 


05:36:24.23 


-06:18:03.6 


21.33 ±0.16 


17.01 ± 0.21 


14.96 ±0.10 




171 


05:36:24.48 


-06:22:23.2 


15.60 ± 0.02 


12.68 ± 0.01 


10.98 ± 0.02 




172 


05:36:24.61 


-06:22:41.7 






15.67 ± 0.04 




173 


05:36:24.63 


-06:20:11.8 


22.77 ± 0.48 




16.09 ± 0.06 




174 


05:36:24.72 


-06:21:54.2 


21.13 ±0.10 


17.05 ±0.18 


15.05 ± 0.01 




175 


05:36:24.83 


-06:25:07.5 


18.91 ± 0.06 


14.72 ± 0.04 


12.65 ± 0.03 




177 


05:36:25.28 


-06:23:07.2 


16.88 ± 0.03 


14.67 ± 0.01 


13.84 ±0.01 




179 


05:36:25.41 


-06:24:31.2 


16.14 ±0.02 


13.90 ±0.01 


13.03 ± 0.01 




181 


05:36:25.85 


-06:23:31.4 




18.08 ± 0.20 


16.09 ± 0.02 




182 


05:36:25.86 


-06:20:55.2 






17.88 ± 0.07 




183 


05:36:25.88 


-06:24:58.6 




16.10 ±0.04 


10.80 ±0.01 
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No. 


RA 


Dec 


I 


/ 








/~>AAA\ 

(2000) 


/TAAA\ 

(2000) 


mag 


mag 


mag 


Comment 


184 


05:36:26.08 


-06:26:15.5 


17.41 ± 0.02 


14.22 ± 0.03 


13.00 ± 0.03 




185 


05:36:26.10 


-06:24:51.7 


16.26 ± 0.04 


14.72 ± 0.01 


13.84 ±0.03 




186 


05:36:26.16 


-06:20:14.8 


22.58 ± 0.41 


17.32 ± 0.22 


14.17 ± 0.06 




187 


05:36:26.71 


-06:26:29.0 


12.65 ± 0.01 


11.04 ±0.03 


10.00 ± 0.02 




188 


05:36:26.83 


-06:24:57.5 


15.05 ± 0.02 


13.23 ±0.01 


12.26 ± 0.01 




189 


05:36:26.93 


-06:18:17.5 




17.04 ±0.18 


15.10 ±0.11 




190 


05:36:27.01 


-06:21:07.6 


16.93 ± 0.02 


15.55 ± 0.07 


14.71 ±0.01 




192 


05:36:27.72 


-06:23:12.3 


14.18 ±0.02 


11.97 ±0.01 


10.83 ± 0.01 




194 


05:36:27.75 


-06:21:57.3 




16.96 ± 0.06 


14.90 ± 0.01 




195 


05:36:27.89 


-06:25:36.0 


14.58 ±0.01 


12.93 ± 0.03 


12.04 ±0.02 




196 


05:36:27.96 


-06:26:44.6 




18.61: 


15.40 ±0.14 




197 


05:36:28.35 


-06:22:30.6 






18.41 ± 0.11 


Ext/Gal 


198 


05:36:28.43 


-06:23:04.8 






19.58 ± 0.45 


Ext/Gal 


199 


05:36:28.66 


-06:24:08.7 






17.82 ± 0.07 


Ext/Gal 


201 


05:36:29.71 


-06:22:52.9 






17.78 ± 0.06 




202 


05:36:29.89 


-06:24:29.2 






16.92 ±0.11 


H 2 knot 


204 


05:36:30.10 


-06:23:10.2 


14.56 ± 0.01 


12.17 ±0.03 


10.28 ± 0.01 




205 


05:36:30.50 


-06:23:56.6 


17.15 ±0.02 


15.07 ± 0.04 


14.13 ±0.01 




207 


05:36:31.26 


-06:23:49.9 






17.65 ± 0.04 




208 


05:36:31.84 


-06:23:23.1 


16.88 ± 0.02 


14.94 ± 0.04 


13.82 ±0.01 




210 


05:36:32.52 


-06:19:21.2 


12.43 ±0.11 


11.03 ±0.03 


9.91 ± 0.03 


Double 


212 


05:36:32.68 


-06:21:47.0 






16.78 ± 0.03 




213 


05:36:32.91 


-06:22:33.7 






16.38 ± 0.03 




214 


05:36:33.06 


-06:21:05.0 






16.06 ± 0.02 




215 


05:36:33.81 


-06:19:27.0 


15.40 ± 0.01 


13.44 ± 0.03 


12.57 ± 0.02 
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Table 11. List of all YSO candidates with intrinsic IR excess. (65 sources) 



No. a 






Intrinsic IR excess 






Spectral features 


Extinction / HR-pos 


YSO method 1 " Remarks 




3.6 


4.5 


5 8 


8.0 


14.3 




Ha Li 


A / (mas) 




23 




X 


X 


X 








1.69 /PMS 


IV 


29 




X 














V 


33 




X 














V 


46 




X 


X 


X 


X 




em 


0.00 /PMS 


II 


56 




X 


X 


X 






em 


0.24 /PMS 


II 


60 








X 








1.23 /MS 


III ZAMS ? 


62 




X 














V 


64 






X 










0.71 /PMS 


IV Close to MS 


66 a 


X 


X 




X 


X 




em 


1.55/ PMS 


II 


66 b 


X 


X 




X 


X 




em 


2.08 / PMS 


II 


68 




X 


X 


X 






em 


0.28 / PMS 


II 


73 




X 


X 


X 


X 




em 


0.97 / PMS 


II 


75 a 








X 


X 




em abs 


1.39 /PMS 


I 


75 b 








X 


X 




em 


1.80 / PMS 


II 


104 




X 


X 


X 










V 


109 






X 


X 










V 


110 a 








X 


X 








V 


110 b 








X 


X 








V 


113 






X 


X 










V 


114 




X 


X 


X 


X 




em abs 


0.18 /PMS 


I 


115 








X 










V 


116 








X 










V 


117 


X 


X 










em 


0.03 / PMS 


ii 


119 




X 


X 








em 


0.20 / PMS 


ii 


121 




X 


X 


X 


X 




em 


1.26 /PMS 


ii 


122 


X 


X 


X 


X 


X 








V 


123 






X 


X 










V 


124 




X 


X 


sat 


X 








YSO with an outflow* 1 


126 






X 


X 










V 


132 






X 


X 










V 


133 




X 


X 


X 


X 






4.66 /PMS 


IV 


134 




X 


X 


X 


X 








V 


137 


X 


X 


X 


X 










V 


138 a 








X 


X 






-4.61 


IV 


138 b 








X 


X 






-3.46 


IV 


140 


X 


X 


X 


X 








1.14/MS 


IV ZAMS? 


141 






X 


X 








-1.40 


IV 


142 


X 


X 


X 


X 


X 




em 


3.84 /PMS 


II 


144 a 


X 


X 




X 


X 


X 


em abs 


0.53 /PMS 


I 


144 b 


X 


X 




X 


X 


X 


em 


0.47 / PMS 


II 


145 


X 


X 


X 


X 


X 




em 


3.28 / PMS 


II 


152 


X 


X 


X 


X 


X 




em abs 


1.85 /PMS 


I 


153 






X 


X 










V 


156 






X 










-0.11 


IV 


157 


X 


X 


X 










-0.97 


IV 


162 






X 


X 


X 








V 


167 








X 








3.11 /PMS 


IV Close to MS 


169 






X 


X 


X 








V 


171 




X 


X 


X 






em 


1.52 /PMS 


II 


172 




sat 


sat 


sat 


X 








YSO with an outflow* 1 


174 


X 


X 


X 


X 








2.07 / PMS 


IV 


175 






X 


X 


X 






2.18 /PMS 


III 


179 






X 


X 






em 


0.13 /PMS 


II 


181 


X 


X 


X 


X 










V 


183 


X 


X 


X 


X 


X 








V 


185 


X 


X 


X 










0.09 /PMS 


III 


188 






X 


X 






em abs 


0.28 / PMS 


I 


192 




X 


X 


X 


X 




em 


0.58 / PMS 


II 


194 


X 


X 


X 


X 










V 


195 








X 


X 




em abs 


0.62 / PMS 


i 


204 


X 


X 


X 


X 


X 




em abs 


0.34 /PMS 


i 


205 






X 


X 






em 


0.24 / PMS 


ii 


208 






X 


X 






em 


0.53 /PMS 


ii 


210 a 




X 




X 


X 




em abs 


0.24 / PMS 


i 


210 b 




X 




X 


X 




em 


0.06 / PMS 


ii 



a Bold source numbers are used to indicate intrinsic IR excess. 
b See Galfalk & Olofsson 120071 . 
Reliability class of YSO candidate: 

I: Ho- in emission and Li I /I6707 in absoiption. 

II: Ha in emission and intrinsic IR excess. 

Ill: Non-zero extinction, stellar spectra and located on a PMS track. 
IV: Non-zero extinction, point-like and located on a PMS track. 
V: Intrinsic IR excess only. 
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Table 12. List of all additional YSO candidates (no intrinsic IR excess). (24 sources) 



No. 


Intrinsic IR excess 


Spectral features 


Extinction / HR-pos 


YSO method" 


Remarks 




3.6 4.5 5.8 8.0 14.3 K s 


Ha 


Li 


A / fmae) 






15 




em 


abs 


0.20 /PMS 


I 




28 








1.37 /PMS 


IV 




30 a 




em 




0.33 / PMS 


III 




31 








2.79 / PMS 


IV 




47 








1.46 /MS 


IV 


ZAMS? 


55 




em 


abs 


0.03 / PMS 


I 




61 




em 


abs 


0.01 / PMS 


I 




65 




em 




1.38 /PMS 


III 




78 




em 


abs 


0.07 / PMS 






80 








3.41 / PMS 


IV 




92 








2.11 / PMS 


IV 




118 








2.28 / PMS 


IV 




120 




em 




0.46 /PMS 


III 


Close to MS 


128 








3.45 / PMS 


IV 




148 








5.05 / PMS 


IV 




160 




em 


abs 


0.68 / PMS 






168 








4.37 / PMS 


IV 




170 








1.87 /PMS 


IV 




177 




em 


abs 


0.06 /PMS 






184 




em 




0.73 / PMS 


III 




186 








4.12 /PMS 


IV 




187 




em 


abs 


0.39 / PMS 


I 




190 








1.16 /PMS 


III 


Close to MS 


215 




em 


abs 


0.09 / PMS 


I 





Reliability class of YSO candidate: 

I: Ha in emission and Li I A6101 in absorption. 

Ill: Non-zero extinction, stellar spectra and located on a PMS track. 

IV: Non-zero extinction, point-like and located on a PMS track. 
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Hu« [erg/s/omZ/«] 




Fig. 24. Optical spectra - Exposure times and / band magnitudes are given for each source. 
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Fig. 25. Evidence of youth - Ha in emission and Li I /16707 in absorption. The shifts in wavelength for each source is an artefact 
caused by each star being differently positioned within the wide slit (1'.'2). 
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